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Abstract The epoch of galaxy assembly from 2 < 2 < 4 marks a critical stage during the 

Cn ' evolution of today's galaxy population. During this period the star-formation activity in the 

> 



Universe was at its peak level, and the structural patterns observed among galaxies in the local 
Universe were not yet in place. A variety of novel techniques have been employed over the 



t?^ ' past decade to assemble multiwavelength observations of galaxies during this important epoch. 

o : 

In this primarily observational review, I present a census of the methods used to find distant 
galaxies and the empirical constraints on their multiwavelength luminosities and colors. I then 
discuss what is known about the stellar content and past histories of star formation in high- 
redshift galaxies; their interstellar contents including dust, gas, and heavy elements; and their 
structural and dynamical properties. I conclude by considering some of the most pressing and 
open questions regarding the physics of high-redshift galaxies, which are to be addressed with 
future facilities. 
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1 INTRODUCTION 

Understanding the detailed formation and evolution of the galaxies we observe to- 
day remains one of the great challenges of modern cosmology. An exceedingly rich 
variety of galaxy properties exists in terms of luminosity, mass, color, structure, 
gas content, heavy-element enrichment, and environme nt, many of which, in turn, 



Blanton fc Moustakas f 20091 ). 



are strongly correlated with each other. As reviewed by[ 
based on the latest generation of wide-field surveys of the local Universe, as- 
tronomers have constructed an exquisitely detailed and statistically robust de- 
scription of the galaxy population today. These results are crucial in terms of 
providing a boundary condition, or endpoint, for our description of the formation 
and evolution of galaxies. Ultimately, we strive to tell that story from beginning 
to end. To do so, we must assemble additional observations and theories. 

Multiple complementary approaches can be used to construct the history of 
galaxy formation. These include ab initio analytic models or numerical simula- 
tions; examinations of the fossil record contained in the ages, metallicities, and 
phase-space distributions of stars in nearby galaxies; and direct observations of 
distant galaxies, for which the cosmologically significant lookback time allows 
a probe of the Universe at an earlier time. In order to test theoretical mod- 
els of galaxy formation at every time step, observations of galaxies over a wide 
range of lookback times are required. Furthermore, a robust translation must be 
performed between the observer's empirical quantities of luminosity, color, and 
velocity dispersion, and the theorist's physical quantities of stellar and dynamical 
mass, current star-formation rate and past star-formation history. 

In the comparison between observations and theoretical models of galaxy evolu- 
tion, an important recent development is the establishment of a precision cosmo- 
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logical framework. Observations of the cosmic microwave background radiation, 

large-scale structure, Type la supernova, the abundance of galaxy clusters, and 
the expansion rate of the Universe, all appear to be well described by a cosmolog- 
ical model in which the Universe is spatially flat, with the dominant component 
of the mass-energy density in the form of dark energy, and the remainder con- 
sisting mostly of cold dark matter with a small fraction of baryons. The initial 
spectrum of density fluctuations in this model is adiabatic, Gaussian, and nearly 
scale invariant. The most recent determination of cos mological paramete rs from 



the Wilk inson Microwave Aniso tropy Probe (WMAP) (jSpergel et al. 20031 ) is pre- 



sented in 



Komatsu et al. (2011. ) , and highlights the fact that most parameters are 



determined with better than 5 — 10% precision. 

Constraining the background cosmological parameters is a crucial part of un- 
derstanding galaxy formation not only for converting apparent quantities such 
as flux and angular size, respectively, into intrinsic ones such as luminosity, 
star-formation rate, stellar mass, and physical size. Cosmological parameters 
are also required for precise theoretical calculations whose predictions can be 
compared with observations, because, according to leading models, the underly- 
ing set of cosmological parameters determine how tiny primordial dark matter 
density fluctuations evolve under the influence of gravity into the large-scale 
spatial distribution of matter in the current Universe. In this framework, the 
collapsed perturbations of dark matter - dark matter halos - serve as the very 
sites of galaxy formation. Based on recent determinations of cosmological pa- 
rameters, massive num erical simulations of the growth of dark matter structure 



have been performed (jBovlan-Kolchin et al. 2003 . 



Springel et al. 2005t ). In or- 



der to compare more directly with observations of galaxy formation, models 
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(either numerical hydrodynamic or semi-analytic) describing the baryonic pro- 
cesses of gas cooling, star formation, and metal enrichment are also required, and 
these, too, are advancing with increased spatial resolution and complexity (e.g.. 



Ceverino. Dekel &: Bournaud 201C 



Somerville et al. 20081 ). 



In order to probe the origin of the global patterns observed in the current 
galaxy population, we must look back to a time before these trends were already 
in place. Furthermore, the old stellar populations of nearby early-type galax- 
i es in dense environm ents suggest that the bulk of their stars formed at z > 2 



(jThomas et al. 20051 ). Therefore, catching the formation of spheroids "in the act" 
requires observations at such early times. Given the strong correlation between 
the properties of the spheroidal components of galaxies and their central black 
holes, the epoch when the spheroids are forming holds special interest for explain- 
ing this connection. As described in more detail during the course of this review, 
the redshift range, 2 < z < 4, corresponding to a lookback time of ~ 10 — 12 Gyr, 
is ideal for directly observing the progenitors of today's fairly luminous spheroidal 
and disk galaxies while in the very process of attaining the properties that come 
to define them over the next 10 Gyr up to the present day. Towards the end of 
this redshift range, the overall level of "activity" in the Universe - both in terms 
of star formation and black hole accretion - was at its peak value. In contrast to 
what is observed in the present-day Universe, a significant fraction of the most 
massive galaxies still sustained active star formation. Furthermore, the Hubble 
Sequence of disk and elliptical galaxies was not yet in place, and the abundance 
of rich clusters was vanishingly small. The early Universe looked drastically dif- 
ferent from its current state. Therefore, studying this epoch can yield important 
clues about the evolution of galaxies. 
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Within the past decade, there has been incredible progress in the study of 

galaxies in this important redshift range of galaxy assembly. A variety of novel 
techniques have been used to identify distant galaxies, and the sample of galaxies 
with spectroscopic redshifts at 2 < z < 4 now numbers well into the thou- 
sands. Although we are far from approaching the overwhelming statistical power 
of the giant local redshift surveys such as the Slo an Digital Sky Survey (SDSS; 



Abazaiian et al. 2009 



vey (2dFGRS; 



Blanton fc Moustakas 20091 ) and 2dF Galaxy Redshift Sur- 



Colless et al. 20011 ). in terms of number of galaxies surveyed spec- 



troscopically, volume probed, and data quality, key features of the galaxy pop- 
ulation at these early times are emerging. The physical properties inferred for 
distant galaxies have provided important inputs and daunting challenges to state 
of the art theoretical models of galaxy formation. As we look forward to the 
next generation of instrumentation on current and future large ground-based 
telescopes, as well as the James Webb Space Telescope (JWST) in space, and 
ever more sophisticated galaxy formation models, it is worth reviewing what is 
known about the physical properties of high-redshift galaxies at 2 < z < 4. 

This primarily observational review is constructed as follows. In Section [2l 
we provide an overview of the many different techniques that have recently been 
employed for identifying high-redshift galaxies. We continue in Section [3] by 
reviewing the global multiwavelength distributions in luminosity and color for 
galaxies in this redshift range. In Section [H we delve into the techniques used 
to transform empirical quantities such as luminosity and color into physical ones 
relating to galaxy stellar populations. In particular, we focus here on stellar con- 
tent and the history of star-formation activity, as well as the relationship between 
these quantities. Sections [5] and [6] in turn consider what is known about the inter- 
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stellar medium (ISM) of distant galaxies - where the interstellar contents include 

gas, dust, and metals - and their structural properties and dynamics. Although 

many gaps remain in our knowledge of these fundamental physical properties, as 

described in Section [71 future facilities and instrumentation will guide us in our 

quest to assemble a comprehensive picture of the galaxy population during this 

distant yet intriguing epoch in the history of galaxy formation. 

2 HIGH-REDSHIFT GALAXY SELECTION TECHNIQUES 

In this section, we provide a brief historical context for the recent dramatic de- 
velopments in the study of high-redshift galaxies, as well as reviewing several of 
the most common and complementary techniques for identifying distant objects. 

2.1 Historical Context 

Over the past 10-15 years, the study of high-redshift galaxies has truly exploded, 
with an increasing number of surveys for systems at lookback times of order 10 
Gyr. The number of galaxies with spectroscopic redshifts at z > 2 is now well into 
the thousands, and the number whose multiwavelength photometric properties 
identify them as such is more than an order of magnitude larger. Due to the dif- 
ficulties of obtaining optical and near-infrared (near-IR) spectra of faint objects, 
photometric redshifts have played an increasingly common role in describing the 
properties of distant galaxies - bringing with them both the advantages of much 
larger samples, and also the drawbacks of larger uncertainties in derived galaxy 
properti es. Also, with some exceptions like the VIRMOS VLT Deep Survey 
(VVDSi lLe Fevre et al. 20051 ). which targets galaxies for spectroscopy down to a 



given optical magnitude limit, these new results for the most part utilize several 
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novel techniques for efficiently identifying distant galaxies with minimal contam- 
ination by systems at lower redshift. Although effective at isolating high-redshift 
galaxies, all of these selection methods suffer from incompleteness with respect 
to a sample defined in terms of physical quantities such as stellar or dynamical 
mass, or star-formation rate. In this section, we describe the landscape of galaxy 
surveys that have contributed thus far to our picture of the high-redshift galaxy 
population. 



2.2 Rest-frame Ultra- Violet Selection 

One method for selecting distant galaxies is based on their rest-frame UV colors. 
This method was first applied at z ~ 3, and specifically exploits the combined 
effects of neutral hydrogen opacity within a star-forming galaxy and along the 
line of sight through the intergalactic medium (IGM). Accordingly, z ~ 3 star- 
forming galaxies with moderate amounts of dust extinction (less than a factor of 
100 in the rest-frame UV) will have distinctive colors in a UGTZ filter system, with 
fairly flat G — TZ and extremely red U — G colors. The "Lyman Break Technique" 
has been used to ident ify thousands of ga l axies (so-called Lyman Break Galaxies, 



or LBGs) at z ~_3 



Steidel et al. 2003 



19961 ). and, using different filter sets, 
at z ~ 4 and 5 (|Ouchi et al. 20041 ). At z ~ 6 and beyon d, even redder sets of 



three filters have been u sed to identify star- forming galaxies (JBouwens et al. 20071 . 



201C 



Oesch et al. 2010l ). In these latter cases, however, the main spectral break 
between the bluest and middle filter arises due to hydrogen Lya opacity in the 
IGM, as opposed to opacity at the Lyman limit. Rest-frame UV selection using 
the initial set of UGTZ filters has also been extended down to lower redshift 



(jAdelberger et al. 20041 . 



Steidel et al. 20041 ). where galaxies at 1.4 < z < 2.5 are 
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isolated due to a lack of significant spectral break. Their fairly flat rest-frame 

-lya forest line blan keting in the 



Steidel et al. (20041 ). provides an 



UV colors, modulated at 2 < z < 2.5 only by 

[/-band, also prove distinctive. Figure [H from 

illustration of the rest-frame UV selection criteria in UGTZ color space, tuned to 

find galaxies at 1.5 < z < 3.5. 

Given the criterion of detection in the rest-frame UV, the techniques described 
above necessarily select galaxies with ongoing star formation and are not sensi- 
tive to passive galaxies. Furthermore, the windows in color-color selection space 
exclude galaxies whose rest-frame UV continuum shape is indicative of significant 
dust reddening. Ground-based rest-frame UV surveys are also typically charac- 
terized by a rest-frame UV (observed optical) flux limit. Detection in rest-frame 
UV (observed optical) bands results in objects well suited to optical spectroscopic 
follow-up. Accordingly, successful spectroscopic follow-up of high-redshift galax- 
ies has been weighted towards rest-frame-UV-selected samples, probing galaxies 
with active ongoing star formation. 

2.3 Rest-frame Optical/Infrared Selection 

Other techniques are tuned to select galaxies on the basis of their rest-frame 
optical colors and are based on a detection at observed near- or mid-infrared 
(mid-IR) wavelengths. For these techniques, the relevant spectral break is either 
the Balmer break at ~ 3650A, which arises when the integrated stellar spectrum 
from a galaxy at these wavelengths indicates the overall spectral shape of A-stars, 
or else the 4000 A break, which reflects the absorption from ionized metals in the 
atmospheres of late- type stars. Although both breaks are indications of maturity 
in stellar populations, they are by no means equivalent in terms of the underlying 
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stellar populations that cause them. The Balmer break appears in stellar popula- 
tions featuring ongoing star- formation over sustained timescales (> 100 Myr), or 
post-starburst populations 0.3 — 1 Gyr since the cessation of star formation. The 
4000 A break is strongest in passive stellar populations in which the current level 
of star formation has been negligible for more than 1 Gyr. Isolating z > 2 galaxies 
on the basis of their rest-frame opt ical breaks requires deep, near- IR photometry. 



As described in 



Franx et al. (;2003l ) and 



van Dokkum et al. (20031 ) , a threshold of 



Jvega — Kyega > 2.3 is effective at identifyin g objects that doin inate the high 



mass regime of the z > 2 stellar mass function (JKriek et al. 2008bl ). Figure [2] (left 



from lFranx et al. 20031 ) demonstrates the sensitivity of the Jvega — Kyega color to 
mature stellar populations at z > 2. These Distant Red Galaxies (DRGs) in fact 
typically have significa nt dust obscuration (Ay > 1) and active star- formation 



rates (> lOOM© yr ^) (IPapqvich et al. 2006), althoug h some show little evidence 



for ongoing star formation ( 



van Dokkum et al. 20081 ). Regardless, the red rest- 



frame UV to optical colors of DRGs down to current near-IR limits indicate 
typical stellar masses in excess of 10^^ Mq. The main limitation of such studies is 
the limited amount of spectroscopic follow up, due to the optical faintness (typ- 
ically R > 25) of DR Gs and the difficult y of obtaining large samples of near-IR 



spectra (but see, e.g.. 



Kriek et al. 2008bl ). 



Another common technique for isolating galax ies at 1.4 < z < 2.5 consists of 



the so-called ''BzK" method (JDaddi et al. 20041 ). As shown in Figure [2] (right) 



BzK galaxies are identified in if -selected samples of galaxies on the basis of their 
colors in the z — K versus B — z plane. Galaxies with fairly blue B — z colors and 
red z — K colors are selected as star- forming z ~ 2 systems (^^sBzK^^), due to the 
presence of a Balmer break. At the same time, quiescent systems {^^pBzK") at the 
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same redshift are identified on the basis of red colors in both B — z and z — K. Al- 
though fairly general in terms of selecting both star-forming and quiescent galax- 
ies, this method misses the youngest star-forming galax ies at z ~ 2, which lack a 



significant Balmer or 4000 A break ( Reddv et al. 20051 ). While the star- forming 



sBzK galaxies have rest-frame UV colors that are redder on average than the 
corresponding C/l/-selected galaxies at z ~ 2, there is significant overlap between 
these two photometric selection technique down to a fixed i^-band magnitude 
limit - much more than b etween the UV-sele cted and DRG samples, which over- 



lap at only the 10% level (JReddv et al. 20051 ). Other surveys tuned to find high 



redshift galaxies based on their rest-fram e optical or near-IR 



the Gemini Deep Deep Survey (GDDS; 



aroperties include 



Abraham et al. 20041') and the Galax y 



Mass Assembly ultra-deep Spectroscopic Survey (GMASS; 



Cimatti et al. 20081 ). 



2.4 Submillimeter / Mid-Infrared Selection 

The overall increase in the level of star- formation activity in the Universe at earlier 
times results in an increased abundance of extreme, bolometrically ultra-luminous 
systems. These systems emit the bulk of their radiation at rest- frame far-IR wave- 
lengths, because of copious amounts of dust obscuring their star-formation and 
AGN activity. Systems with bolometric luminosities greater than L = IQ^'^Lq 
are commonly referred to as Ultra-luminous Infrared Galaxies (ULIRGs). Sub- 
millimeter and mid-IR instrumentation sensitive to cool and hot dust, respec- 
tively, have enabled the identification of these extreme high-redshift ULIRGs 
on the basis of their reprocessed emission, while multiwavelength imaging and 
spectroscopic follow-up has elucidated the range of their properties. The Sub- 
millimetre Common-user Bolometer Array (SCUBA) has been used to identify 
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Small. Ivison fc Blain 19971 ). 



bolometrically-luminous submllllmeter galaxies (SMGs; 

and the largest set of spectroscoplcally-confirmed such objects was obtained by 

following up the subset of sources with both Fsso/^m > 5 m Jy and Very Large Ar- 



ray ( VLA) 1.4 GHz fluxes greater than F-^ ^ GHz ~ 30/iJy (jChapman et al. 20031 . 



20051 ). The radio fluxes were used to obtain the precise (1" —2") positions required 
for spectroscopic observations, which were not achievable with the coarse (15") 
SCUBA beam. These luminous SMGs with radio counterparts have a median 
redshift of z = 2.2, although, given that the requirement of a radio flux detection 
recovers ~ 50% of the Fgso^m > 5 mJy population, their redshift distribution 
may not be fully representative of the luminous SMG population as a whole. 

The Multiband Imaging Photometer for Spitzer (MIPS) onboard the Spitzer 
Space Telescope has been used to identify high-redshlft dusty so urces on the basis 
of t heir brightness at 24Atm and faintness at optical wavelengths. 



Yan et al. f2007l ) 



and lDev et al. (20081 ) present such samples, identified using slightly different cri- 
teria, but both based on similar criteria of detection at 24;um, with a large ratio 
of mld-IR to optical flux. The majority of these objects are at 1.5 < z < 3, with 
comparable space densities (~ 10~^ Mpc^ ) to those of SMGs. On the other 
hand, the 24/xm-selected sources are characterized by warmer dust temperatures 
and a higher frequency of AGN signatures at mld-IR wavelengths than SMGs 
selected at longer wavelengths (see Section 15. 2p . Quantifying the relative contri- 
butions of star-formation and AGN activity in powering the extreme luminosities 
of both mld-IR and submlllimeter-selected ULIRGs will enable us to isolate the 
underlying nature of these sources. Another key goal consists of constraining the 
relative Importance of major mergers and smooth mass accretion as triggers for 
the ULIRG phase. 
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2.5 Narrowband Selection 

In contrast to the identification of liigh-redshift objects on the basis of broad- 
band spectral shape, the use of a narrowband filter tuned to the redshifted 
wavelength of a specific emission line is effective at isolating objects with large 
emission-line equivalent widths. The most common emission line for which nar- 
rowband filters are designed is hydrogen Lya. As shown in Figure [3] (from 



Gronwall et al. 20071 ). based on images through both the narrowband filter and a 



broadband filter close in wavelength, objects with red broadband minus narrow- 
band colors are flagged as Lya line emitters (LAEs). More th an 2000 LAEs have 



been identifled with ground-based facilities at 2 < z < 8 (e. g. 



Gronwall et al. 20071 . 



Nilsson et al. 2011 



Quchi et al. 2008 



Cowie fc Hu 19981 . 



Rhoads et al. 200d ). 



The star-fornnng LAEs tend to be significantly fainter on average than the UV- 
continuum-selected objects described in Section 12.21 and therefore offer a probe 
of the faint end of the luminosity function. On the other hand, the faint nature of 
these objects leads to a challenge in assembling high signal-to-noise (S/N) mul- 
tiwavelength imaging and spectra for individual LAEs (with information other 
than a measurement of Lya emission) , hindering the determination of their rela- 
tionship to other galaxy populations at similar redshifts. 

3 EMPIRICAL PROBES OF THE HIGH-REDSHIFT GALAXY 
POPULATION 

Before reviewing what is known about the stellar and interstellar content of dis- 
tant galaxies, we must consider the empirical measurements from which these 
physical properties are inferred. The observables here are distributions in lumi- 
nosity and color, which offer some of the most basic and fundamental probes 
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of a galaxy population. In addition to traditional optical photometry trac- 
ing the rest-frame ultraviolet, our view of the global photometric properties 
of galaxies at z > 2 is now based on deep, near-IR surveys from the ground 
using wide-field imagers on 4-meter-class telescopes (e.g., KPNO/NEWFIRM, 
UKIRT/WFCAM, CTIO/ISPI, Palomar/WIRC) and narrow- and wide-field im- 
agers on 8 - 10-meter class telescopes (e.g., Keck/NIRC, VLT/ISAAC, Sub- 
aru/MOIRCS, VLT/HAWK-I), and mid- and far-IR surveys using Spitzer and 
Herschel in space. Multiwavelength observations have, therefore, granted us a 
window into the luminosity and color distributions of high-redshift galaxies span- 
ning from the rest-frame UV through the rest-frame far-IR. 

3.1 Luminosity Functions 

The galaxy luminosity function offers constraints on the overall abundance of ob- 
jects, as well as the integrated luminosity density at a given wavelength. As such, 
the luminosity function provides a key observational baseline for each redshift at 
which it is measured. 

3.1.1 REST-FRAME UV LUMINOSITY FUNCTIONS Some of the 
first luminosity functions to be measured for high-redshift galaxies w ere based on 



optical observations of rest-frame UV-selected z ~ 3 and z ~ 4 LBGs (jSteidel et al. 1993 ) 
probing rest-frame wavelengths of A ~ 1700 A. This work highlighted the im- 
portance of dust extinction for converting the rest-frame UV luminosities of 
star-forming galaxies into unobscured values, in order to obtain dust-corrected 
star- formation rates. We will return to this point in Sections H] and 15. H when 
we consider the star-formation rates and dust content of high-redshift galaxies. 
Rest-frame UV luminosity functions for large samples of distant galaxies have now 
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been estimated by se veral different groups, from z ~ 2 all the way out to z ~ 8 



(JBouwens et al. 20ld ). As in the local Universe, the galaxy luminosity function 



using the Sche chter form of a power 



at high redshift is typically parameterizec _ 

law multiplying an exponential function (jSchechter 19761 ) . The associated free 
parameters to constrain are the characteristic luminosity, L* (or M* in the space 
of absolute magnitude), the f aint-end slope, a, an d th e overall normalization, $*. 



At z ~ 2 — 3, the work of 



Reddv et al. (200Si ) and 



Reddv k Steidel (20091 ) is 



based on the largest set of spectroscopic redshifts in the literature, and, including 
a consideration of the systematic variation of d ust reddening with UV luminos- 
ity, is accordingly the most robust. Specifically, iReddv &: Steidel (2009. ) utilizes 
> 2000 spectroscopic redshifts, and ~ 31,000 2; ~ 2 — 3 photometric candidates 
in 31 independent fields over 0.9 deg . In this paper, the standard LBG selec- 
tion limit of 7^ = 25.5 was extended to fainter magnitudes (~ O.IL*) to obtain 



tighter constraints on the faint-enc 
ruinations at z ~ 4 include those by 



slope, a. Recent lum inosity function deter- 



Bouwens et al. (20071 ). using 4671 S-dropout 



galaxies selected over 580 arcmin^ with deep Hubble Space Telescope {HST) Ad- 
va nced Camera for Surveys (ACS) imaging, down to Af^yAB ~ ~16, and 



by Ivan der Burg. Hildebrandt &: Erben (20101 ). based on ~ 36,000 ^r— dropout 
galaxies selected in the CFHT Legacy Survey in four independent 1 deg^ fields 
down to MjjyjyQ = —18.7. Both of these z ~ 4 surveys, however, are solely based 
on photometric selection, without spectroscopic confirmation. Figure H] reviews 
recent determinations of the rest-frame UV luminosity function at 2 < 2; < 4. 

The rest-frame UV luminosity functions of star-forming galaxies at z ~ 2 — 4 
are characterized by steep faint-end slopes (a ~ —1.6 1.7) and characteristic 



luminosities of M 



AB 



-21. These steep faint-end slopes are in contrast to 



16 Shapley 

the flatter one {a = —1.22) determined from the local far-U V (FUV) luminosit y 



function using Galaxy Evolution Explorer (GALEX) data (jWvder et al. 20051 ). 
Also, the characteristic luminosities at 2; ~ 2 — 4 are roughly three magni- 
tudes brig hter than those determ ined from the local GALEX FUV luminosity 



function. iReddy &: Steidel (20091 ) review other recent luminosity function de- 
terminations at z ~ 2 — 3, highlighting some of the discrepancies in the liter- 
ature, including those determinations with signi ficantly fiatter faint-end slopes 



(jGabasch et al. 2004 . ISawicki fc Thompson 200(1 ). and/or larger counts at the 



brightest luminosities ( Paltani et al. 20071 ) . To address some of these differences 



Reddy Sz Steidel (2003 ) point out the pitfalls associated with attempting to mea- 
sure the luminosity function over small areas (< 50 arcmin ), making incorrect 
assumptions about the nature of the intrinsic mean and dispersion in colors of the 
faintest galaxies, and insufficiently accounting for contamination by low-redshift 
interlopers. Again, we emphasize the importance of spectroscopy for understand- 
ing the redshift selection function of the objects for which the luminosity function 
is being constructed, and the proper characterization of systematic effects mod- 
ulating the volume probed by a given galaxy survey. 

3.1.2 REST-FRAME OPTICAL AND NEAR-INFRARED LUMINOS- 
ITY FUNCTIONS While the rest-frame UV luminosity of distant galaxies 
reflects the emission from massive stars, longer wavelengths probe different as- 
pects of galaxy stellar populations and dust content. Specifically, the rest-frame 
optical luminosity function is more reflective of older stars, although the extent 
to which emission at these wavelengths reflects the integra ted stellar mass de- 



pend s in detail on the star-formation history of the galaxy (jShaplev et al. 2001 



2OO5I ) . The rest-frame optical (T^-band) luminosity function was flrst determined 
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at z ~ 3 bv lShaplev et al. (2001. ). based on the LBG 7^-band luminosity function 
and the distribution oi TZ — Kg colors for a sample of 118 LBGs, 81 of which 
had spectroscopic redshifts. These measurements yielded a steep faint-end slope 
of a = —1.85, which is significantly different from the local determinations with 
a ~ —1, and a characteristic l uminosity of M^ = — 22.98, which is 1.5 magnitudes 



brighter than the local value (JBlanton et al. 20031 ) 



Most recently, based on a much larger set of ~ 1000 i^-band measurements at 
2 < z < 3.5, selected o yer a total area of 378 a rcmin^ from multiple near-IR sur- 
veys of varying depths. iMarchesini et al. (20071 ) construc ted B,V , and -R-band lu - 



Marchesini et al. (20071 ) 



minosity functions. Using mainly photometric redshifts 
find that the faint-end slopes of these luminosity functions are, within the errors, 
consistent with the fairly shallow slope determined for the local optical lumi- 
nosity function, whereas the char acteristic magnitudes are significantly brighter 



(> 1 magnitude). In contrast to 



Shaplev et al. r200lh . 



Marchesini et al. (20071 ) 



find values of a ranging from —1. 1.4. Howev er, it is worth pointing out 

that, in the region of overlap, the 



Marchesini et al. 



ity function for "blue" {J — K < 2.3) galaxies and the iShaplev et al. 



rest-frame F-band luminos- 
rest-frame 



y-band LBG luminosity function are entirely consistent. Both total and "blue" 



y-band luminosity functions 
the luminosity function from 



:"rom 



Marchesini et al. (20071 ) are shown along with 



Shaplev et al. (200ll ) in Fi gure [5l reproduced from 



Marchesini et al. (20071 ) 



Marchesini et al. (20071 ). Potential limitations of the 
analysis are in the small fraction of spectroscopic redshifts (~ 4%) and corre- 
sponding reliance on photometric redshifts, and the small area over which the 
photometry is deep enough to robustly probe the faint end (~ 25 arcmin ). 
With ultra-deep i^-band [Kumit ~ 23 (Vega)] surveys over significantly larger 
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areas (~ 1000s of arcmin^), such as the UKIDSS Ultra-Deep Survey (UDS), the 

much-needed robust constraints on the faint-end slope of the rest-frame optical 

luminosity function will be within reach. Extensive spectroscopic follow-up is 

also necessary t o avoid some of the bias es related to photometr ic redshifts that 



are described in 



Reddv fc Steidel ^20091 ^ 



Reddv et al. (2008l ) and 
In principle, the rest-frame near-IR luminosity is even more closely tied to stel- 
lar mass than the rest-frame optical (despite some of the uncertainties we will dis- 
cuss in Section 14. 3p . The largest study to date of the evolving rest-frame K-hand 
luminosity function is based on the UKIDSS UDS Fir st Data release, presen ted by 



Cirasuolo et al. (201Q ). updating their earlier work (jCirasuolo et al. 2007i ). This 
study features a K-l-z-band-selected catalog of ~ 50, 000 galaxies over 0.7 deg (~ 
10,000 of which are at z > 1.5), which is complete down to i^ = 23 (AB). Given 
the location of the UKIDSS UDS in the Subaru/XMM-Newton Deep Survey 
field, this survey also benefits from extensive multiwaveler igth coverage spanning 



from FUV (GALEX) to mid-IR (Spitzer) wavelengths. 



Cirasuolo et al. (20ld ) 



estimate photometric redshifts and rest-frame i^-band luminosities based on 



spectral-energy distribution (S 
their sources. The depth of the 



D) fits to the multi wavelength photometry of 



Cirasuolo et al. (2010l ) sample is not sufficient to 



trace the faint-end slope and its evolution past 2; ~ 1, so the evolution to high 
redshift is simply quantified in terms of the luminos ity function normaliza tion 



Cirasuolo et al. (20ld ) re- 



and characteristic luminosity. From z ~ to z ~ 2, 
port a brightening in M* by ~ 1 magnitude and a decrease by a factor of ~ 3.5 
in normalization, $*. Again, deeper near-IR photometry and extensive spectro- 
scopic follow-up for near-IR-selected catalogs will be required to probe the full 
evolution of the rest-frame near-IR luminosity function at high redshift. 
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3.1.3 REST-FRAME MID-INFRARED AND BOLOMETRIC LUMI- 
NOSITY FUNCTIONS At even longer rest-frame wavelengths, we begin to 
probe the direct emission from dust. Building on the work of earlier missions 
such as the Infrared Astronomical Satellite (IRAS) and the Infrared Space Ob- 
servatory (ISO), which traced the evolution of infrared-luminous (IR-luminous) 
sources out to z ~ 1, Spitzer has played a crucial role in tracing the global dust 
emission properties for large samples of z > 1 star-forming galaxies. The most 
widely- used tool in this endeavor is the 24;um channel of the MIPS instrument, 
which probes a rest-frame wavelength of A ~ 8/im at z ~ 2, sensitive to the 
emission from polycyclic aromatic hydrocarbon (PAH) emission. Ground-based 
submillimeter observatories such as SCUBA have also played a key role in under- 
standing the evolution of most IR-luminous sources, but for much smaller galaxy 
samples. Overall, out to z ~ 1, the evolution of the mid-IR and total IR luminos- 
i ty function is charact erized by both significant luminosity and density evolution 



(JLe Floc'h et al. 20051 ). IR-luminous galaxies were more numerous in the past, 
and the IR luminosity density at z ~ 1 is dominated by luminous infrared galaxies 
(LIRGs), with Lr p > lO^^Lc^. 



Recent work by 



Perez- Gonzalez et al. (20051 ) 



Caputi et al. f2007l ^. and 



Rodighiero et al. (20ld ) 



have characterized the luminosity f unctions of MIPS 24 jUm- selected galaxies at 



z > 1. We focus here on the results of 



Caputi et al. f;2007l ') and 



Rodighiero et al. f20ld ) 



who used more conservative criteria for excluding AGNs from their samples, and 
adopted conversions between rest-frame 8/um and total IR luminosities that are 
most consistent with observed constraints (and agree well with each other). These 
studies select galaxies down to a flux limit of S'(24;um) = 80/iJy in the 0.08 deg of 
the Great Observatories Origins Deep Survey (GOODS) North and South fields. 
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while iRodighiero et al. (2010. ) additionally include a shallower catalog down to 



a limit of S'(24^m) 
shown in Figure [U 



400// Jy in the .85 deg^ of the VVDS-SWIRE area. As 



Caputi et al. (2007i ) determine mid- and tot al IR luminos 



Rodighiero et al. (2010l ) construct the 



ity functions at z ~ 1 and z ~ 2, while 

corresponding luminosity functions in nine redshift bins from z ~ to z ~ 2.5. 

Both of these works echo the previously determined increase of an order of 
magnitude in the bolometric IR luminosity density from z ^ to z ^ 1, ari d the 
heightened abundance of both LIRGs and ULIRGs. IRodighiero et al. (2010i ) find 
that the IR lumino sity density is roughly constant from z ~ 1 to z ~ 2.5, while 



Caputi et al. (2007i ) find evidence for a slight decline. For both of these studies, 



it is worth pointing out that their completeness limits at z ~ 2 in 8^m lumi- 
nosity translate into bolometric IR luminosities of ~ IO^^Lq. Therefore, these 
luminosity functions only directly probe the ULIRG regime. Furthermore, nei- 
ther study can place constraints on the faint-end slope of the luminosity function, 
and instead adopt a fixed parameter of a = —1.2 based on local observations. 
Using indirect estimates of the rest-frame 8/um and bolometric IR luminosity func- 
tions at z -^ 2, inferre d from the extinctio n-corrected luminosities of UV-selected 
star-forming galaxies, iReddv et al. (20081 ) extend the determination of the IR lu- 



minosity function into the LIRG regime, and suggest that a significantly steeper 
faint-end slope may be required. Accordingly, ULIRGs make a sub-dominant 
(~ 25%) contribution to the far-IR luminosity density at z ~ 2. Future direct 
observations of this fainter IR regime with Herschel, and even more sensitive 
planned facilities such as the Single Aperture Far-Infrared (SAFIR) observatory, 
will be crucial in untangling the evolution of the IR luminosity density and the 
history of obscured star formation in the Universe. 
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3.2 Color-Magnitude Diagrams 

A standard tool for describing the galaxy population in the local Universe is the 
optical "color-magnitude" diagram, in which rest-frame optical colors are plotted 
as a function of rest-frame optical absolute magnitude. The added dimension 
of color proves very useful for separating different types of galaxies. Indeed, 
while the division of the galaxy population into different types (e.g., late- type 
spirals and irregulars, and early- type ellipticals and lenticulars) , and the corre- 
lation between galaxy structure and color have been well-known for a long time, 
the advent of large spectroscopic surveys such as SDSS have allowed for a quan- 



titative description of the bimoda 
incredibly robust statistics. Both 



ity in galaxy photom etric properties, based on 



Strateva et al. (20011 ) and 



Baldrv et al. (20041 ) 



have presented striking evidence that local galaxies occupy a bimodal distribu- 
tion in the space of u — r colors, and that the fraction of bluer galaxies increases 
at fainter optical absolute magnitudes. Galaxies occupying the so-called "red- 
sequence" part of the bimodal distribution primarily consist of morphologically 
early-type galaxies and follow a very tight relationship between color and mag- 
nitude in which more luminous galaxies have redder colors. At the same time, 
objects residing in the "blue cloud" are predominantly late-type galaxies, and 
show a looser, though still systematic, variation of color and luminosity in the 
same sense. The bimodal distribution of galaxy colors as a function of luminos- 
ity, along with the correlation of color and other galaxy properties, suggests two 
distinct types of formation histories. It is therefore of critical interest to trace 
the galaxy color-magnitude diagram to earlier times, and to determine how far 
back the bimodality in the galaxy distribution persis ts. 



Based on the COMBO-17 survey. 



Bell et al. (2004 ) demonstrate that bimodal- 
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ity is detected in the color-magnitude diagram out to 2; ~ 1. Now, using the re- 
sults from recent near- and mid-IR-selected surveys, other groups have considered 
the question of galaxy bimodality at even higher redshifts. As shown in Figure [71 
with a sample of 1021 Spitzer/IRAC 4.5^m-selected objects from the GMASS 
survey with optical thr ough mid-IR SEDs, 190 of which have spectroscopic red- 



shifts above z = 1.4, 



Cassata et al. (20081 ) demon strate that a bimoda. 



Brammer et al. (2009 ) uses 



ity in 



galaxy rest-frame U — B colors persists up to z = 2. 
^ 25,000 objects with K < 22.8 (AB) selected from the NEWFIRM Medium- 
Band Survey (NMBS), and detect a bimodality in rest-frame U — V colors out to 
z ~ 2.5. The increasing importance of obscured star formation at higher redshifts 
ten ds to cause contaminat ion of the red sequence by dusty, star-forming galaxies, 



and iBrammer et al. (20091 ) show a much cleaner division in dust-corrected colors 
or when considering only galaxies whose mid-IR flux limits of S'(24^m) < 20/iJy 
indicate the presence of little or no dust. Based on a smal ler sample of 28 K - 



selected galaxies with spectroscopic redshifts 2 < z < 3, iKriek et al. (2008a 



detect a red sequence in rest-frame U — B colors, quantified in terms of a signif- 
icant over density of galaxies within a narrow bin of red rest-frame U — B color. 
These high-redshift red-sequence galaxies are characterized by little or no ongo- 
ing star formation, with strong Balmer breaks indicating that they are likely in 



a post-sta rburst phase. Whl 



sample in 



e suggestive of the existence of a red sequence, the 



Kriek et al. (2008al ) is too small to test for the presence of bimodality 



in gala xy colors. 



Both 



Kriek et al. f2008al ') and 



Cassata et al. (20081 ) stress the need for precise 



(and ideally spectroscopic) redshift information when attempting to detect fea- 
tures in the galaxy rest-frame color distribution such as bimodality or the presence 
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of a red sequence. While they lack actual spect roscopic confirmation, galaxies in 



the NMBS presented by 



Brammer et al. (2003 ) have extremely accurate photo- 



metric redshifts, with Az/{l + z) < 0.02 at z > 1.7. Photometric redshifts typical 
of other studies, with errors oi Az/{l + z) ~ 0.1 or worse, will lead to random and 
systematic errors in inferred rest-frame UV and optical colors that will wash out 
these trends when a single color is considered. However, using only photometric 



Williams et al. f2009l ) 



redshifts but considering two rest-frame colors together, 
discern a bimodal behavior in the space of rest-frame U — V versus V — J color- 
color space out to z ~ 2.5. The rest-frame colors for this study are inferred from 
optical, near-IR, and Spitzer/IRAC photometry drawn from the UKIDSS UDS, 
the Subaru- XMM Deep Survey (SXDS), and the 5piteer Wide- Area Infrared Ex- 
tragalactic Survey (SWIRE), and effectively separate quiescent, non-star- forming 
galaxies from their actively star-forming counterparts. The evolving locations in 
color space and very existence of the "red sequence" and "blue cloud" provide 
important constraints on models of galaxy formation. 

4 THE STAR-FORMATION RATES AND STELLAR CON- 
TENT OF HIGH-REDSHIFT GALAXIES 

While the luminosities and colors of galaxies represent basic and fundamental 
observables, we seek to translate these measurements into physical quantities. 
Specifically, the SEDs of galaxies are commonly interpreted in terms of the cur- 
rent rate of star formation and its past history, as well as the integrated stellar 
content of galaxies. In this section, we summarize both the simple methods used 
to infer star-formation rates from specific luminosities, as well as the techniques 
used to model the stellar populations of high-redshift galaxies based on mul- 
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tiwavelength SEDs. Along the way, we highlight the results of applying these 

methods to determine the global history of star formation, the build-up of stellar 

mass density, and the relationship between star-formation rate and stellar mass 

in distant galaxies. 

4.1 Star-formation Rate Indicators 

As reviewed by 



Kennicutt (19981 ). there are several diagnostics of star- formation 
activity in external galaxies. In the study of distant galaxies at z > 2, integrated 
light measurements tracing young stellar populations are used to infer the rate at 
which stars are being produced. These measurements include rest-frame UV and 
IR luminosities, hydrogen recombination emission-line luminosities, and stacked 
X-ray and radio luminosities (due to the sensitivity limits of current X-ray and 
radio facilities, z > 2 star-forming galaxies are typically only detected in a sta- 
tistical sense). In general, these diagnostics are only sensitive to the presence of 
massive stars. Therefore, an estimate of the total star-formation rate requires 
the assumption of a particular form for the stellar initial mass function (IMF), 
which is then used to extrapolate down to the low stellar masses that dominate 
the integrated stellar mass. 

The most commonly used star-formation diagnostic for high-redshift galaxies 
is the ultraviolet continuum luminosity over the wavelength range 1500 — 2800 A, 
which is dominated by O and B stars and directly related to the star-formation 
rate in galaxies where star formation has been proceeding at a roughly constant 
rate on ~ 10®-year timescales. The rest-frame UV luminosity at 2 < z < 4 is 
based on observed-frame optical photometry. As described in Section 15. H the 
colors of star-forming galaxies at high redshift suggest that dust extinction leads 
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to a significant attenuation of tlie rest-frame UV luminosity. Therefore, the rest- 
frame UV luminosities based on optical apparent fluxes must be corrected for 
dust when calculating intrinsic star-formation rates. 

Dust absorbing the radiation from massive stars is heated and re-radiates this 
absorbed energy in the far-IR region of the spectrum. Therefore, the far-IR lu- 
minosity from dust also provides a tracer of the rate of massive star form ation. 



The range "far-IR" is commonly defined as 8 — 1000/im (JKennicutt 19981 ). and 
the luminosity spanning this wavelength range is adopted as a tracer of star 
formation for systems in which star formation has proceeded continuously for 
at least 10^ — 10^-year timescales. For galaxies with bolometric IR luminosities 
large enough {Ljr > IO^^Lq) to be detected with ground-based submillimeter 
telescopes such as SCUBA, the observed 850/xm flux can be converted into a 
far-IR luminosity with an assumption of dust temperature (see Section r5.2p . For 



z ^ 2 galaxies with smaller bolometric luminosities (down to Lj^ ~ 10 Lg 



0^ 



deep Spitzer/MIPS 24/um observations have provided a proxy for far-IR lumi- 
nosity. At z ~ 2, the 24^m channel probes a rest wavele ngth of S/xm, the loca - 



tion of strong PAH em i ssion. Local galaxy templ ates from 



Charv fc Elbaz f200lh 



Elbaz et al. f 20021 ). and 



Dale fc Helou f2002l ) 



Rieke et al. (20091 ) are used to 



convert from the mid-IR to total IR luminosity, which corresponds to an increase 
by an order of magnitude. Recent Herschel Photodetector Array Camera and 
Spectrometer (PACS) 160//m measurements of the far-IR luminosities of star- 
forming galaxies at z ~ 2 suggest that MIPS 24/im p roxies for far-IR lum inosi- 



ties tend to yield overestimates by factors of ~ 4 — 7.5 (jNordon et al. 2010l ). This 
apparent discrepancy will require additional study. Spitzer/MIPS 70/im data 
probing rest-frame 20 — 25/im at z ~ 2 have not been used as widely because of 
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the limit of the depth of the existing Spitzer data, in which very few individual 

z > 2 galaxies are detected. 

As photons from massive stars ionize nearby interstellar gas, the emission from 
recombining hydrogen gas in these star-forming regions serves as a proxy for the 
rate of production of ionizing photons, and, by extension, the formation rate of 
massive stars. The ionizing flux is dominated by emission from the most mas- 
sive (> IOMq) stars and provides an estimate of the instantaneous rate of star 
formation. With ground-based observations. Ha emission lines have been mea- 
sured for star- forming galaxies up to z ^ 2.6 (at which redshift the thermal 
backg round begins to dominate the noise) and used to estimate star formation 



rates (JErb et al. 2006bl . 



Forster Schreiber et al. 2009 ). The weaker H/3 feature 



has also been used to estimate star-forma tion rates at redshifts beyond z = 2.6 

). In addition to Balmer lines, the Lya 



(jMannucci et al. 20091 . 



Pettini et al. 2001 



feature has been used as a proxy for star-formation rate. However, due to the 
resonant nature of the Lya transition, this line is especially sensitive to the ef- 
fects of dust extinction and scattering, which can both lead to the preferential 
destruction of Lya photons and their diffusion over a large area with reduced 
surface brightness. Accordingly, there is evidence that even dust-corrected obser- 
vations of Lya from individual galaxies tend to underpredict the star formation 
rate, when co mpared with othe r indicators such as rest-frame UV luminosity or 



Ha emission ( Haves et al. 2011 



Additional star-formation rate diagnostics are based on X-ray and radio lumi- 
nosities. High-mass X-ray binaries, young supernova remnants, ai id hot interstel- 

lar gas contribute to the X-ray luminosity in star-forming galaxies. iRanalli. Comastri &: Sett 
calibrate the relation between rest-frame 2 — 10 keV and far-IR luminosities for 
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such systems, leading to a relation between hard X-ray luminosity and star- 

jetween far-IR and 1.4 GHz 



Yun. Reddv fc Condon f200lh 



formation rate. Based on the tight linear correlation 
radio luminosity among nearby star-forming galaxies, 
derive the relation between star-formation rate and radio luminosity. The deep- 
est current Chandra 0.25 — 2.0 keV (observed) and VLA 1.4 GHz radio data (e.g., 
in the GOODS-N region) are still not sufficient for detecting all but the most 
luminous star-forming galaxies and AGNs at z > 2. Therefore, most studies 
that make use o f X-ray or radio estimates of star-forniation rates use stackiri g 



techniques (e.g.. 



Daddi et al. 2007I . 



Pannella et al. 20091 . 



Reddv k Steidel 20041 ) 



yielding only sample-averaged properties. 



4.2 Evolution of the Star- formation Rate Density 

The evolution of the star- formation-rate density is considered one of the most fun- 
damental observational descriptions of the galaxy population as a whole. Match- 
ing and explaining this observed evolution is often used as a benchmark of success 
for theoretical models of galaxy formation. The evolution of the star-formation- 
rate density is constructed by integrating the luminosity function at a specific 
wavelength sensitive to star formation (e.g., rest-frame UV, Hq, far-IR, or radio), 
in order to obtain the associated luminosity density. Then, a conversion between 
luminosity and star-formation rate (as described in Section 14. ip is used to ob- 
tain the associated star-formation-rate density. In such studies, it is crucial to 
specify th e limits down to which the luminosity function is integrated at each red- 
shift (e.g.. lBouwens et al. 20071 ). Also, at shorter wavelengths, such as rest-frame 



UV and Ha, corrections for dust extinction must be appli ed in order to obtai n 
the unobscured star-formation rate density (see Section [5. UlCalzetti et al. 2000l ). 
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These corrections are sometimes applied in an average sense to the integrated 

luminosity density, or else in a luminosity-dependent fashio n, taking into a ccount 



the relationship between luminosity and dust obscuration (JHopkins 2004 ). 



The firs t attempts to chart the st ar-formation histo ry of the Universe were pre 



sented by 



Madau et al. (19961 ) and 



Lilly et al. (19961) . At th e pre sent, an exten 



Hopkins f2004l 'l and 



Hopkins &: Beacom (200( 



sive set of measurements has been compiled by[ 
from z ~ to z ~ 6, using a variety of star-formation rate indicators includ- 
ing rest-frame UV, far-IR, radio, and X-ray continuum luminosities, and both 
Balmer and forbidden [Oil] emission lines. These measurements were all con- 
verted to a common cosmology {Qm = 0.3, Qa = 0.7, Hq = 70 km s^ Mpc^^), 



an d rest-frame UV and optica 
In 



Bouwens et al. (20071 . 



measurements were corrected for dust extinction. 
2010l ). rest-frame UV luminosity functions, both cor- 
rected and uncorrected for dust extinction, are used to map the evolution of the 
star-formation-rate density to z ~ 8. As shown in Figure [51 in both of these 
compilations, the evolution of the star-formation-rate density is characterized by 
an order -of-magnitude i ncrea se from z ^ to z ^ 2. In more detail, the compila- 



tions bv iHopkins (20041 ) and 



Hopkins fc Beacom (2006: 1 show the dust-corrected 



star-formation-rate density increasing by an order of magnitude from its local 
value already by z ^ 1 and remaining roughly flat to z ~ 2. The evolution in 



Bouwens et al. (20071 . 



20 id ) is characterized by a smooth rise all the way from 
2; ~ to z '~ 2, due to a different extinction correction adopted for the rest- 
frame UV data at z < 2. In both versions of the cosmic star-formation history, 



the global star-formation-rate density remains rough 
and z ~ 4, and then declines towards higher redshift. 



y constant between z ^ 2 



Hopkins &: Beacom (20061 ) 



propose that the star-formation history is constrained to within '^ 30 — 50% up 
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to z ~ 1 a nd within a factor of ~ 3 at higher redshifts. Additionally, recent 

results from 



Reddv et al. (20081 ) suggest that ~ 70 — 80% of the star-formation- 
rate density at z ~ 2 — 3 is produced by galaxies with bolometric luminosities 
Lboi < lO^^-^0, and that the highly-o bscured ULIRGs selected by submillime- 



ter surveys (e.g., 



Chapman et al. 20051 ). although individually luminous, do not 



dominate the star-formation-rate density. 

4.3 Stellar Population Synthesis Models 

While luminosities tied to individual specific wavelength ranges are commonly 
used to infer the current rate of star formation, the multiwavelength SED of a 
galaxy can be interpreted in terms of its integrated stellar, dust and metal con- 
tent, and its past history of star formation. This technique is referred to a s stellar 

1! 
Tinslev (1968i ). In the 

intervening decades, stellar population synthesis models have become increasingly 
sophisticated, and have played a crucial role in inferring the physical properties of 
galaxies both near and far based on their photometric and spectroscopic proper- 
ties. Population synthesis models are also used to estimate photometric redshifts, 
in so far as the observed multiwavelength photometry of a galaxy is fit in terms 
of not only best-fitting stellar population parameters, but also redshift. 

The basic ingredients of stellar population synthesis models consist of a the- 
oretical prescription for all stages of stellar evolution as a function of mass and 
metallicity, as well as a stellar spectral library of the observed properties of stars 
at different positions in the Hertzsprung-Russell diagram. The time evolution of 
the integrated spectrum from a coeval population of stars described by a specific 
stellar IMF and metallicity can then be predicted. For a given star-formation 
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history, the evolution of the integrated spectrum of stars is computed by sum- 
ming the contributions of simple stellar populations formed at successive time 
steps, with the normalization at each time step determined by the correspond- 
ing star-formation rate. Thus, stellar population synthesis models can predict 
the evolution of the integrated spectrum of a stellar population for an arbitrary 
star- formation history, IMF, and stellar metallicity. Stellar population s ynthesis 



models are often coupled with theoretical mo dels (jCharlot &: Fall 200d ) or em- 
pirical parameterizations (jCalzetti et al. 2000i ) of dust extinction to explain the 
observed properties of galaxies in terms of both stellar and dust content. The re- 
sulting model galaxy spectrum can be passed through photometric filters in order 
to predict the evolution of luminosities, colors, and mass-to-light ratios in specific 
photometric bands. Additional properties of the stellar population can also be 
calculated, including the rate of Type la and Type II supernovae explosions, the 
mass in stellar remnants and recycled gas, and the ionizing photon luminosity. 

Currently, there are several different stellar population synthesis models em- 
ployed for fitti ng the SEDs of distant g alaxies. Perhaps most widely used are 



the models of 



Bruzual fc Chariot (2003. ) and their current version (Chariot &: 



Bruzual 2011, in preparation), though other significant theoretical efforts in 



elude the PEGASE and Starburst99 models (JFioc fc Rocca-Volmerange 19971 . 



Leitherer et al. 19991). and, ni ore recently, the models of 



Maraston (20051 ) and 



Conrov. Gunn &: White (20091 ). There is much ongoing discussion in the litera- 
ture regarding the systematic uncertainties and differences among these different 
stellar population synthesis codes, in terms of their descriptions of various stages 
of stellar evolution off of the main sequence. 

There has been a particular focus on the treatment of the thermally-pulsating 
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asymptotic giant branch (TP-AGB). TP-AGB stars are red giants with low- to 

intermediate-mass main sequence progenitors, which make a significant contri- 
bution to the near-IR luminosity of a simple stellar p opulation at ages between 



0.5 and 2.0 Gyr. As emphasized by 



Maraston (20051 ). derived properties such 



as the age and near-IR mass-to-light ratio, and, correspondingly, the inferred 
stellar mass, will be very sensitive to the treatment of TP-AGB stars in systems 
in an evolutionary state that prominently features these stars. As a result, the 
stellar populations in which the correct treatment of the TP-AGB matters the 
most are post-starburst systems, as o pposed to those in which star formation is 
proceeding at a roughly constant rate (|Daddi et al. 20071 ) . As shown in Figure [H 



in the models of 



Maraston (20051 ). TP-AGB stars make a much more significant 



contr ibution to the integrated 



els of 



galaxy luminosity at ~ 1 Gyr than in the mod- 



Bruzual &: Chariot (20031 ). resulting in systematical 



masses and ages, in particular for post-starburst systems (JMaraston et al. 20061 ). 



y lower derived stella r 



Rather than adopting a specific prescr iption for the contribution of the TP-AGB 



phase, IConrov. Gunn &: White (20091 ) parameterize the luminosities and effec- 
tive temperatures of TP-AGB stars as variables to be constrained by the actual 
data. Systematic uncertainties in TP-AGB parameters therefore translate into 
systematic errors in other der ived properties suc h as stellar mass and past star- 



formation history. Recently, 



Kriek et al. (2010l ) constructed a composite SED 



over the rest-frame range 1200 — 40000 A for a sample of 62 post-starburst galax- 
ies at 0.7 < z < 2.0 d rawn from the NEWF IRM Medium-Band Survey. In 



contrast to the results oflMaraston et al. (20061 ). this entire SED is fit well by the 



Bruzual &: Chariot (2003i ) models, whereas the rest-frame optical and near- 



gions of the spectrum are not simultaneously fit by the models of 



Maraston 



Rre- 
. The 
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SED-fitting results of 



3 



r 



Kriek et al. (20101 ) suggest that the iMarastoni models give 
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too much weight to the TP-AGB phase. Clearly, consensus has yet to be reached 
about the influence of TP-AGB stars in the integrated spectra of galaxies. 

Uncertainties in the detailed nature of the stellar IMF also translate into sys- 
tematic uncertainties in the conversion between luminosity and stellar mass. For 
the range of stellar populations typically observed at high redshift (i.e., when 
the Universe was less than a few Gyr old), and over the rest-frame UV to 
near-IR wavelength range where current observations prob e, models with the 



same star-form ation history and the assumption of either a 



Chabrier f2003l ') or 



Salpeter (1955l ) IMF over the mass range 0.1 — lOOMm pr oduce virtually identical 



colors. However, the model assuming a 



Chabrier (20031 ) IMF corresponds to a 



stellar mass a factor of ~ 1.8 lo wer. Although bo th IMFs are described by power- 



Chabrier (20031 ) IMF follows a log-normal distri- 



law functions at M > IMq, the 

but ion below IMp) , turning over at the so-called "characteristic mass," whereas 



the 



Salpeter (19551 ) IMF continues to increase as a power-law all the way down 
to O.IMq. 



I 

Evidence in the local and low-redshift Universe suggests that the .Chabrier (20031 ) 



IMF yields stellar M/L ratios for lower-mass elliptical galaxies (with velocity 
dispersions of a ~ 200 km s~ ) in agreement with their luminosities and dynam- 
ical mass estimates. On the other hand, elliptical galaxies with larger dynamical 



masses appear to be characterize d by either steeper (i.e. , more Salpeter-li 



or higher dark-matter fractions ([Graves &: Faber 2010l . 



ke) IMFs 



Treu et al. 20ld ). There 



are also extreme environments such as the Galactic Center and surrounding star- 
burst clusters suc h as the Arches, i n which direct evidence for IMF variations 



has been reported (jStolte et al. 20021 ). At higher redshift, based on more indirect 
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methods such as measuring the simultaneous U—V color and stellar i?-band M/L 

ratio evolution for elliptical galaxies atz~0to2;~0.8 ( with the latter estim ated 



van Dokkum (20081 ) in 



from the observed evolution in the fundamental plane), 
fers that the stellar IMF for elliptical galaxy progenitors had a flatter slope at 
z ~ 4 in the regime near IMq. Furthermore, the characteristic mass at which 
the log-normal portion of the IMF turns over is inferred to shift to ~ 2Mq. 
Adopting a completely independent approach based on a comparison of the star- 
formation r ates and stellar masses of vigorously star- forming galaxies at z ~ 2, 



Dave (20081 ) also infers a higher characteristic mass at this redshift. An upward 
shift in char acteristic mass r esults in an IMF that is more "bottom light" than the 
present-day IChabrier (20031 ) functio n, and a lower conversion fac tor from light 



to stellar mass. In contrast, recently 



van Dokkum &: Conrov (20ld ) analyzed the 



spectra of four elliptical galaxies in the Virgo cluster and, based on the strength 
of stellar absorption features tracing M < O.3M0 stars, concluded that the IMF 
in the massive star-forming progenitors of these systems had a "bottom heavy" 
I MF steeper than Salpeter betwee n 0.1 and I.OMq. 

Bastian. Covev &: Mever (201Q ) review evidence for IMF variations both in 



the local Universe and at significant cosmological distances for the most part 
with skepticism. Clearly additional study is required to quantify or rule out the 
evidence for IMF variations as a function of environment. Stellar population 
synthesis models featuring explicit parameterizations of the most unconstrained 
phases of stellar evolution, as well as the form of the stellar IMF, will yield con- 
fidence intervals for derived stellar population parameters that more accurately 
reflect the systematic uncerta inties associated with stellar population modeling 



(IConrov. Gunn fc White 20091 ). 



34 Shapley 

In spite of the significant uncertainties associated with stellar population syn- 
thesis modeling, it is now standard practice to use such models to infer basic 
physical properties of galaxies over a wide range of redshifts. In the local Uni- 
verse for a sample of > 10^ SDSS galaxies, a combination of optical broadband 
photometric properties and spectral indices (the 4000 A spectral break and the 



strength of Ba 
thesis code of 



mer absorption lines) hav e been modeled using the population syn- 



Bruzual fc Chariot (20031 ) to infer galax y properties such as star - 



formation history, dust attenuation, and stellar mass (JKauffmann et al. 2003al ). 
At higher redshifts, stellar population modeling is typically tuned to broadband 
photometry alone, as rest-frame optical spectra are not of sufficient quality to 
measure stellar absorption features or detailed continuum shape at high S/N. At 
z > 2, stellar population synthesis modeling is only possible if the rest-frame 
SED is probed at wavelengths both above and below age-sensitive spectral dis- 
continuities such as the Balmer or 4000 A break. Therefore, optical photometry 
probing the rest-frame UV must be, at the very least, combined with data probing 
the rest-frame optical regime (i.e. near-IR observed wavelengths), and, prefer- 
ably, with additional Spitzer/IRAC photometry probing t he rest-frame near-IR . 
Early models of z > 2 stell ar populations were featured in 



Shaolevet al. f200l!). and 



Sawicki fc Yee f 19981 ) 



Papovich. Dickinson &: Ferguson (2001. ) (LBGs) and 



Forster Schreiber et al. f 20041) fPRGs). A fair ly recent example of z -^ 2 stellar 



population modeling from 



Muzzin et al. f2009l ) is shown in Figure [TOj This anal- 



ysis is based on extremely well-sampled SEDs for spectroscopically-confirmed 
galaxies including optical, near-IR, and Spitzer/IRAC photometry as well as 
binned Gemini/GNIRS near-I R spectra, and f e atures a systematic compa rison 



of the parameters derived from 



Maraston_(2005|), 



Bruzual fc Chariot ^20031 ). and 
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updated Chariot & Bruzual models. Stellar population synthesis models are now 

routinely used to derive high-redshift galaxy physical properties. Such model- 
ing is a standard component of the derivation of global distributions such as the 
galaxy stellar mass function (Section 14. 5p . 

For the modeling of high-redshift galaxy stellar populations, in addition to 
the adoption of a particular stellar population synthesis code, the stellar IMF 
and metallicity are assumed parameters. A specific, wavelength-dependent dust 
extinction law (e.g., starburst, SMC, Milky Way) is also assumed. The star- 
formation history, SFR{t), is commonly parameterized in the form of an expo- 
nential decline, with SFR{t) = SFRq x exp(— t/r). In this case, r, the e- folding 
time, and t, the time since the onset of star formation, are both parameters to 
constrain. With such a parameterization, a continuous star-formation history 
corresponds to r = cxd and SFR{t) = SFRq. In addition to single episodes de- 
scribing the star-formation rate as a function of time, more complex functions for 
the history of star formation have been considered. In particular, two-component 
models consisting of the linear combination of an old, high mass-to-light ratio 
component and a younger population with ongoing star formation, have been 
used to constrain how much stellar mass from the old stellar population could be 
"hiding" under the glare of a younger burst of star formation. In addition to r 
and t, the parameters commonly derived for high-redshift galaxies are indicators 
of the degree of dust extinction [E{B — V) or Ay], the current star- formation 
rate, and the integrated stellar mass. For a given stellar population synthesis 
code, the stellar mass has been demonstrated to be the best-constrained param- 

i 1 I 1 I 1 

eter (jPapovich. Dickinson &: Ferguson 20011 : IShaplev et al. 20011 . i2005i ) , whereas 
other parameters are more subject to uncertainties in the nature of the star- 
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formation history (i.e., r), which is difficu 
multiwavelength information. Recently, 



t to constrain in the absence of external 



Maraston et al. (20101 ) have in fact ar- 



gued that so-called "inverted-r" models [i.e. models in which the star-formation 
rate increases with time as SFR{t) = SFRoxexp^+t/r)] provide a better descrip- 
tion of the extinction and star-formation rate based on rest-frame UV data alone, 
as well as the star- formation rates and stellar masses of m ock high-redshift galax - 



Papovich et al. (201 ih 



ies constructed from semi-analytic models. In addition, 
demonstrate that the globally averaged relations between the star-formation rates 
and stellar masses of galaxies at high redshift (see Section 14. 6p appear to favor 
rising star-formation histories, as opposed to ones that are constant or declining. 
The best parameterization of the star-formation histories of high-redshift galaxies 
is clearly still a matter of debate. 

4.4 The Diversity of High-Redshift Stellar Populations 

The methods described above have been used to investigate the stellar popu- 
lations of high-redshift galaxies selected using the various techniques discussed 
in Section [21 Based on the results of stellar population syntehsis modeling, we 
can make some general comments about the range of stellar populations ob- 
served, while keeping in mind the uncertainties inherent to the modeling process, 
and the biases that different selection techniques impose. If the current sam- 
ples of UV-selected galaxies are roughly characterized by star-f ormation rates of 



10 - 100M(T) yr~\ typical stellar masses o f 1 - 5 x lO^^M© (jErb et al. 2006cl . 



Reddv fc Steidel 2004 



Shaplev et al. 2001 



and moderate amounts of extinc- 



tion in the rest-frame UV (factor of ~ 5), the very complementary sample of 
rest-frame optically selected DRGs are typically characterized by higher star- 
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formation rates (> WOMq yr~^), stellar masses (down to the typical survey lim- 
its of A' ~ 21 Vega) of ~ W^^Mq, and larger amo unts of dust extinction in the 



rest-frame UV and o ptical ranges of the spectrum ([Forster Schreiber et al. 20041 . 



Papovich et al. 20061 ) . At t he same time, a minority of DRGs show little evidence 



for ongoing star formation (jvan Dokkum et al. 20081 ) . In addition to their prodi 



gious star-formation rates, SMGs appear to be characterized by stellar masses 
that are comparable to those of the typic al DRGs, and several times larger on av- 



erage than those of UV-selected galaxies (|Borvs et al. 2005i : iMichalowski. Hiorth fc Watson 20ld ) 
However, it is worth keeping in mind that the effects of dust extinction and AGN 
contamination on the broadband SED lead to larger systematic uncertainties in 
the derived stellar populations of SMGs. At the other end of the spectrum, so to 
speak, the stellar populations of LAEs are characteri zed by less dust extin ction on 



average even than those of the UV-selected galaxies (jGawiser et al. 20071 ) . Due to 
the typically faint rest-frame UV luminosities of emission-line selected galaxies, 
the LAEs also tend to be faint in the rest-frame UV continuum, with lower star- 
formation rates than those of UV-selected gala xies, which were targeted down to 



a brighter continuum limit (,Kornei et al. 2010l ). 



In addition to considering the typical properties of galaxies as a function of se- 
lection method - which may not have anything other than historical value - it is 
also worth mentioning the range of star-formation histories observed as a function 
of mass. In particular, at high stellar masses (M > 10^^ Mq, assuming a Salpeter 
IMF from 0.1 — IOOMq), there exist at z ~ 2 not only active star-forming galaxies 
with SFR > IOOM0 yr~^, but also passive, evolved galaxies with little evidence 
for ongoing star formation (jKriek et al. 2008ai . ivan Dokkum et al. 2008 ) . Quies- 
cence and mature stellar populations constitute the physical interpretation of the 
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Kriek et al. (2008a ) and described 



empirically-derived red sequence reported by 
in Section 13. 2i These quiescent galaxies appear to consitute ~ 40 — 50% of the 
most massive galaxies {Mgtar > 10^^ Mq at z ~ 2 — 3). Based on these results, it 
is worth noting that, in contrast to the patterns observed in the local Universe, 
at least half of most massive galaxies at z ^ 2 a re still actively in the process of 



forming fpaddi et al. 2007 



Papovich et al. 20061 ). and that there is an incredible 



diversity observed among the star-formation histories of these massive galaxies. 
At the same time, it is a challenge to explain massive galaxies at early times 
with little evidence for ongoing star formation, given that theoretical models pre- 
dict copious rates of mass a ccretion at z ~ 2 for such massive systems (e.g.. 



Dekel. Sari fc Ceverino 20091 ). 



4.5 Evolution of the Stellar Mass Density 

Stellar population synthesis modeling provides a powerful tool for estimating the 
global evolution of the stellar content in galaxies, which reflects the combined 
processes of star formation and mergers. This evolution is described by con- 
structing the galaxy stellar mass function at a range of redshifts. Analogous to 
the galaxy luminosity function, the stellar mass function is parameterized in terms 
of a characteristic mass, M*^^^., low-mass slope, a, and overall normalization, <I>*. 
At each redshift, the stellar mass function can be integrated to determine the 
corresponding stellar mass density. Alternatively, the growth in stellar mass as 
a function of galaxy mass can provide important insights into galaxy formation 
models. The cosmic stellar mass density should also reflect the integral of past 
star formation in the Universe, and therefore the integral of the star-formation- 
rate density described in Section 14.21 can be compared for consistency with the 
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independently-derived stellar-mass density. 

In the local Universe, the galaxy stellar mass function has been determined 

from large samples of galaxies drawn from the 2dFGRS matched to the 2 Mi- 



cron All-Sky Su rvey (2MASS), and SDSS (JBaldrv. Glazebrook fc Driver 20081 : 



Cole et al. 2001 



Presently, stellar mass functions have been measured out to 
z ~ 5, suggesting that roughly half of the local stellar mass density appears to 
be in place at z ~ 1. Measurements of the stellar mass function at z > 1 are 
based on samples with multiwavelength (optical and IR) photometry and primar- 
ily photometric redshifts. T he first determination o f the stellar mass function at 
z ~ 2 — 3 was presented in [Dickinson et al. (2003l ) , based on a sample of rest- 
frame -B-band-selected objects in the Hubble Deep Field North (HDF-N) with 
both optical and near-IR photometry. Subsequently, many other groups have 
measured stellar mass function s at z > 2, selecti i ig galaxies at optical, near-IR 



and mid-IR wavelengths (e.g. 



Fontana et al. 2004 



Both 



20061: 



Drorv et al. 20051 : 



Pozzetti et al. 2007'). 



Perez- Gonzalez et al. (20081 ) and 



Eisner. Feulner fc Hodd 20081 : 



Marchesini et al. (20091 ) construct stel 



lar mass functions at z >2 based on fairly deep (near-IR and mid-IR magnitude 
limits of 23 — 25 AB) and wide-area (500 — 700 arcmin^) surveys, selected, re- 
spectively, with IRAC and iT-band data. The larger areas of these surveys, 
compared to previous determina tions, provide results th at are less susceptible 



to cosmic variance. As shown in 



Marchesini et al. f 20091 ). the 2 < z < 3 stellar 



mass functions in the literature, when integrated over a fixed range in stellar 
mass (10^ < M.s tor /M(7) < 10^^), vary b y a factor of ~ 3 in stellar mass density. 
Figure [11] (from iMarchesini et al. 2009 ) features a compilation of global stellar 
mass density estimates as a function of redshift. At z ~ 2, the reported fraction 
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of the local stellar mass density that is in plac e ranges from ^ 8 — 25 %, while 



Marchesini et al. (2003 1 also of- 



that number drops to ~ 4 — 12% at z ~ 3.5. 
fer an in-depth analysis of the random and systematic uncertainties involved in 
constructing a stellar mass function from the modeling of multiwavelength pho- 
tometry. These include the errors associated with photometric redshifts, cosmic 
variance , differences among stellar pop ulation synthesis codes (as parameterized 



by, e.g.. IConrov. Gunn &: White 20091 ). choice of stellar IMF, stellar metallicity, 
and extinction law. A proper accounting for these sources of random and system- 
atic error significantly increases the uncertainties on derived quantities such as 
the evolution of the global stellar mass density, as well as the evolution of galaxies 
as a function of st ellar mass. For example , when only taking into account ran- 



dom uncertainties, 



Marchesini et al. (2009 ) find that the abundance of galaxies 



below the characteristic mass evolves more strongly with redshift than that of 
the most massive galaxies {Mstar > lO^"*^'^), which show a lack of strong evolution 
in number density. However, when the full random and systematic error budget 
is accounted for, more significant evolution in the abundance of massive galaxies 
cannot be ruled out. 

Additional uncertainties in the stellar mass density at high redshift result from 
uncertainties in the low-mass slope of the stellar mass function. While most works 
adopt a in the range to —1.0 to —1.4, the stellar mass regime crucial for constrain- 
ing this parameter is n ot well probed with cur rent observations. For example. 



the dataset featured in 

low ~ IO^^^Mq. According to 



Marchesini et al. f2009ll suff ers from incompleteness be- 



Reddv &: Steidel (2003 ) , if the stellar mass function 



has a steeper low-mass slope, as suggested by the steep faint-end slope of the rest- 
frame UV luminosity function and the relationship between UV luminosity and 
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stellar mass, up to ^ 50% of the stellar mass density may be contained in galaxies 

with stellar masses < 10^^ Mp, as opposed to the ^ 1 0—20% inferred from extrap- 



olating the Schechter fits of iMarchesini et al. (20091 ). In fact, using a -ff -selected 



sample based on significantly deeper ne ar-IR imaging with Subaru/MOIRCS 



{K = 24.1 Vega), 



Kaiisawa et al. (2003 ) estimate a steeper low-mass slope of 



a 



-1.5 at z ~ 2 and a = —1.6 at z ~ 3, as well as perha ps a trend of a 



Fontana et al. (20061 ). This result is 



steepening with redshift, also suggested by 
intriguing, but the area over which the low-mass slope is adequately probed is only 
28 arcmin . Clearly, data of this depth must be collected over a significantly wider 
area to minimize the effects of cosmic variance and obtain more robust constraints 
on the stellar mass density at z > 2. In particular, constraining the abundance 
and dust-extinction properties of faint, low-mass galaxies at high redshift will 
prove very important for comparisons of the past integr al of the star-formation - 



rate density with the stellar mass density at each redshift (JReddv fc Steidel 2003 ) . 



Careful comparisons of this sort are crucial for determining whether or not the in- 



tegral of globa 



past star formation indicates a disc repancy with the global stellar 



mass density (jWilkins. Trentham &: Hopkins 20081 ). 



4.6 SFR - Mstar Scaling Relations 

In addition to considering the evolution of star-formation rates and stellar masses 
in galaxies separately, the evolution of the relationship between these quanti- 
ties provides important clues as to how stellar mass builds up in galaxies as 

a function of reds hift and mass. At z -^ 2, this relationship is highlighted by 

I 

Daddi et al. (2007l ) in a sample of star-forming sBzK galaxies (passive pBzK 
galaxies and galaxies with no MIPS 24^m detection were excluded from this anal- 
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ysis). Using star- formation rates estimated from either dust- corrected UV lumi - 



Daddi et al. f2007l ) 



nosity, or the sum of mid-IR and uncorrected UV luminosities, 
find a strong correlation between star-formation rate (SFR) and stellar mass, 
described by the relation, SFR oc M^i^^., and shown in Figure [T2j The ultra- 
luminous SMGs are outliers to the SFR— M star trend, with st ar-formation rates a 
factor of ~ 10 higher than expected, given their stellar masses. iPannella et al. (20091 ) 



used 1.4 GHz radio stacking observations to estimate average star- formation 
rates for non-AGN sBzK objects in the Cosmic Evolution Survey (COSMOS) 
field, binned by stellar mass . The radio stack ing analysis reveals the trend 
SFR oc M^i^^, consistent with .Daddi et al. (2007. ). Accordingly, the specific star- 
formation rate, i.e. the star-formation rate divided by stellar mass, appears to be 
roughly constant over an order of magnitude in stellar mass (10^*^ — IO^^Mq, as- 
suming a Salpeter IMF from 0.1 — IOOMq). Similar trends between star- formation 
rate and stellar ma ss ha ve been observed arn ong star-forming galaxies at z ~ 1 by 



Elbaz et al. (2007i ) and lNoeske et al. (20071 ). but with a lower overall normaliza- 



tion, such that, at a given fixed stellar mass, the expected star-formation rate is 
a factor of '^ 4 lower. A comparison with th e trend observed in the local Universe 



(JBrinchmann et al. 20041. 



Elbaz et al. 2007 ) indicates an evolution by a factor of 

-n — — 



30 - 40 (|Daddi et al. 2007i . iPannella et al. 2009l l. 



On the other hand, perhaps the tightness of the SFR — Mgtar correlation has 
been overemphasized, as a result of selection effects. Using the sam e ultra-deep. 



Kajisawa et al. (2010l ) investigate 



-ff-selected sample described in Section 14.51 

the SFR — Mstar relation from 0.5 < z < 3.5. At z ~ 2, these authors find 

significantly more scatter at the high-stellar-mass end of the relation. In this 



study, UV dust-corrected and mid-IR-|-UV star-formation rates are estimated 
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using the same techniques as in iDaddi et al. (2007i ) , but there is no requirement 

for galaxies to be detected at 24/im. At M^tar = 10^"^ M©, for example, UV dust- 
corrected star-formation rates range from 1 — IOOOMq yr~^. Furthermore, the 
power-law slope between star-formation rate and mass, and therefore, the specific 
star-formation rate, tend to decrease at Mstar > IO^'^'^'^Mq. Given that the slope 
and small scatter of the SFR — Mstar relation has been interpreted in terms 
of lending support to galaxy formation models in which smoo th gas accretion 
dominates the growth of galaxies at high redshift f|Dave 2008i ). it is crucial to 
characterize these quantities accurately and in an unbiased manner. The above 
studies are largely based on photometric redshifts, which tend to increase the 
uncertainties in all derived physical properties. Therefore, larger spectroscopic 
samples are needed. Based on such spectroscopic studies, a description of the 
distribution of star-formation rates as a function of stellar mass for a stellar- 
mass selected sample will provide the ideal observational probe of this potentially 
meaningful trend. 

5 THE INTERSTELLAR CONTENT OF HIGH-REDSHIFT GALAX- 
IES 

The stellar content of galaxies offers an incomplete version of the story of their 
formation and evolution, which must be filled in by a characterization of their 
interstellar environments. Indeed, the multiwavelength study of the current and 
past history of star formation in galaxies cannot be constructed without un- 
derstanding the nature of dust. In addition to attenuating and reddening the 
radiation from stars, dust also reradiates the absorbed emission in the IR with a 
characteristic overall spectral shape that depends on temperature. The detailed 
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emission spectrum in the mid-IR offers a further probe of the energy sources 

heating dust grains (i.e., radiation from star formation or an AGN). While dust 

reprocesses the Ught from stars, it is the cool gas content of galaxies that forms 

the very fuel for star formation. The elevated rate of star formation in galaxies 

in the early Universe is a direct result of the large cool gas fractions in these 

systems. As stars evolve and die, they return gas and heavy elements to the 

ISM. The patterns of chemical enrichment in both gas and stars therefore reflect 

the past history of star formation, and of gas inflows and outflows in galaxies. 

In this section, we review what is known about the interstellar contents of high- 

redshift galaxies. We begin by discussing both dust extinction and reprocessed 

emission, and then turn to observations of the molecular gas and metals, all of 

which provide important insights into the nature of distant galaxies. 

5.1 Dust Extinction 

Starburst galaxies in the local Universe follow a correlation between attenuation 
and reddening. As these galaxies become more attenuated in the rest-frame UV 
(as probed by the ratio of far-IR to UV luminosities), their rest-frame UV con- 
tinua become redder. Based on multiwavelength observations spanning from the 
UV to far-IR, this correlation has been quantified as a relationship between the 
rest-frame UV slope, /3 (where f \ oc A^), and ^ifinn (the attenuation at 1600 A), 



such that ^1600 = 4.43 + 1.99/3 
burst obscuration curve of 



Meurer. Heckman &: Calzetti 19991 ). The star 



1 
Calzetti et al. (200Q ) describes the wavelength depen- 



dence of effective attenuation when dust is distribu ted in a patchy fore ground 



screen, relative to young stars. For reference, the 



Calzetti et al. (2000 ) curve 



lacks the 2175 A bump characteristic of the Milky Way extinction curve, and has ; 
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different ratio of total to selective extinction, with Ry = Ay /E{B — V) = 4.05, a s 



opp osed to the average M ilky Way value of 3.1 (jCardelli. Clayton fc Mathis igsg) 



The 



Calzetti et al. f200Q l law is also "grayer" than the SMC law (jPrevot et al. 19841 1 



rising less steeply in the near-UV, and therefore implying more attenuation for 
a given reddening in the rest-frame UV. This "grayness" likely stems from a ge- 
ometrical configuration between gas and stars in which the dust is distributed 
i n a patchy fo reground-like screen, and some stars suffer little or no extinction 



([Calzetti 200 Ij ). This starburst attenu ation curve predicts an almost identi cal re- 



lation between ^leoo and /3 to that of 



Meurer. Heckman &: Calzetti (19991 ) . with 



the assumption of fairly uniform intrinsic rest-frame UV colors for starburst galax- 



ies. 



Soon after the discovery of LBGs at z ~ 3, their observed r ange of rest-frame 



UV colors was interpreted in terms o f a range of dust red dening ([Calzetti et al. 200d : 



Meurer. Heckman &: Calzetti 19991 : 



Steidel et al. 19991 '). in analogy with the de- 



scription of local starbursts. The ICalzetti et al.l obscuration curve was used to 
translate between rest-frame UV color and E{B — V), and, by extension, ylieoO) 
with the average value of E{ B — V) = 0.15 cor responding to an attenuation 



factor of ~ 4.7 at 1500 A. In 



Steidel et al. (19991 ). this factor was used to cor- 



rect the observed UV-luminosity densities at z ~ 3 and z ^ A and show that 
only a small fraction of the intrinsic UV radiation typically escapes from even 
UV-selected galaxies at high redshift. Further more, consistent with the tren d ob- 



Adelbereer fc Steidel (200d ) find 



served among starbursts in the local Universe, 

that objects with greater bolometric luminosities suffer more extinction in the 

rest-frame UV. However, the trend at z ~ 3 is offset from the local one in the 



sense that, for a given bolometric luminosity (i.e. star- formation rate), the ex- 
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tinction is significantly smaller at higher redshift. iBouwens et al. (2009l ) have 

investigated the average reddening as a function of rest-frame UV luminosity 

(uncorrected for dust), demonstrating that, at z ~ 2.5 and z ~ 4.0, /3 becomes 

bluer for fainter objects. Furthermore, the average /3 value of the most UV- 

luminous galaxies is redder at z ~ 2.5 than at z ~ 4 .0. The trend between and 



UV 
and 



uminosity was not apparent in the datasets of lAdelberger &: Steidel (200d ) 



Reddv et al. (20081 ). most likely due to the smaller dynamic range of UV 



luminosities probed. 
At this point, the 



Calzetti et al. (200d ) law is the most commonly adopted 



extinction law for modeling the stellar populations of high-redshift galaxies (Sec- 
tion 2]). Therefore, it is important to consider the empirical support for this 
choice, based on comparisons between UV-extinction-corrected and extinction- 
free estimates of star- formation rates. While individual star- forming galaxies 
are not detected in deep Chandra X-ray and and VLA 1.4 GHz radio imaging 
observations, stacking methods have proven very pow erful for comparing dif - 



Reddv fc Steidel f200^ 



ferent multiwavelength star- formation rate indicators, 
measure the stacked X-ray and radio fluxes for spectroscopically-confirmed, UV- 
selected star-forming galaxies at 1.5 < z < 3.0 in the GOODS-N field. On 
average, the inferred X-ray and radio-derived star-formation rates are, respec- 
tively, 42 and 56 Mq yr~^. The UV star- formation rate, extinc tion corrected 



Calzetti et al. 



based on rest-frame UV colors and assuming the validity of the 
law, is 50 Mq yr^^ - consistent with the radio and X-ray estimates. Further- 
more, the ratio of star-formation rates derived from X-ray and uncorrected UV 
fiuxes implies a factor of ~ 4.5 — 5.0 attenu ation in the UV, cons istent with 



the attenuation inferred from the UV colors. 



Pannella et al. (20091 ) perform a 
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similar radio stacking analysis for star-forming sBzK galaxies with photometric 

redshifts at z ~ 2 in the COSMOS field, finding that the UV attenuation, ^1500, 

derived from the ratio between the extinction-free radio-derived star-formation 

rate and the UV-uncorrected star- formation rate, agrees well with that inferred 

on the basis of rest-frame UV colors alone, assuming the 

law applies. 

S'piteer/MIPS 24|im obser vations have allowed for a test of extinction on a 



Calzetti et al. f200d ) 



per object basis at z ~ 2. jReddv et al. (2010 



20061 ) compute mid-IR lumi- 



nosities from observed 24;um fluxes for spectroscopically-confirmed UV-selected 
galaxies, and extr apolate these to total IR luminosities using local templates 



(|Elbaz et al. 20021 ). The relation between extinction {Lfir/Liqqq) and redden- 
ing (/3) for these UV-selected galaxies is compared with the relation among local 
starbursts, revealing that the majorit y of objects (thoug h not all - see below) 



follow the local trend. Fur thermore , 



Reddv et al. (20061 ) confirm the result of 



Adelberger &: Steidel (20001 ) of the positive correlation between bolometric lu- 



minosity and extinction, but now based on a Spitzer/MIVQ estimate of the IR 
luminosity. The evolution in extinction with redshift is also reproduced using 
Spitzer, which suggests that extinction in the rest-frame UV is a factor of ~ 10 
smaller at ^ ~ 2 than at 2; ~ 0, at fixed bolometric luminosity. iDaddi et al. (20071 ) 



perform an analogous test of UV extinction laws, using star-formation rates for 
star-for ming sBzK galaxi es derived from both res t-frame UV an d Spitzer 24//m 



fluxes. 



Daddi et al. (20071 ) also conclude that the 



Calzetti et al. 



law is valid for 



the the majority of systems at z ^ 2. 



1 

While the success of the ICalzetti et al. (200Q ) law is impressive, in terms of 

correcting UV luminosities and matching various extinction-free tracers of star- 
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formation rates for large samples of high-redshift galaxies, it is also important 

to highlight the cases in which it appears to fail. For example, in ultraluminous 



SMGs at 2 ~ 2, the dust-corrected UV luminosities (assuming the iCalzetti et al. 



law) underpredict the bolometric luminosities suggested 



ay submillimeter and 



radio luminosities by as much as a f actor of ~ 10 — 100 (jChapman et al. 2005 



Daddi et al. 2007 



Reddy et al. 2006i). A related yet ra ther non-intuitive fact is 



that > 50% of the SMGs in the 



Chapman et al. (20051 ) sai nple have rest-frame 



UV c olors that actually satisfy the UV-selection criteria of iSteidel et al. (20031 . 



20041 ) . The mismatch between predicted (based on rest-UV color) and observed 
IR luminosities, shown in Figure [13] (left), may arise in systems where regions 
of massive star formation are completely opaque to UV radiation, and the UV 
radiation that does escape is from regions that are di sj oint from the dusty ones 



dominati ng the bolometric ou tput (JDaddi et al. 20071 ). 



that the 



Reddv et al. 1^20061 1 find 



Calzetti et al. (2000l ) law appears to break down at z ^ 2 for system s 



more luminous than ~ 2 x IO^^Lq. On the other hand, iMagdis et al. (2010a 



demonstrate that, at z ~ 3, the starburst attenuation curve yields consistent 
results for g alaxies with bolometr ic luminosities as large as IO^^Lq. In general. 



however, the 



Calzetti et al. (2000l ) law does not appear valid for describing dust 



extinction in the most luminous sources at z > 2. 

While the local starburst relation appears to underpredict the UV attenua- 
tion for the most bolometrically luminous systems, there is a discrepancy in the 
opposite sense for objects that, based on their stellar population modeling, ap- 
pear "young," with best-fit ages (assuming constant star-formation histories) less 
than 100 Myr. This second discrepancy was highlighted most dramatically in 
the case of the strongly gravitationally-lensed objects MS1512-cB58 (or "cB58") 



Galaxies at 



■^-r4 



(|Pettini et al. 200Q ) and the Cosmic Eye (jSmail et al. 20071 ) . Both of these ob- 
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jects have rest-frame UV colors and continuum slopes that suggest dust-corrected 
star-formation rates several times larger than what is actually measured using 



Spitzer/MIPS mid-IR photometry and InfraRed Spectrograph ( 



(a fa ctor of ~ 3 — 5 for cB58, and ~ 8 for the Cosmic Eye) (|Siana et al. 20091 . 



RS) spectroscopy 



20081 ). This discrepancy, shown in Figure [T3l (right), was pr eviously noted for 



cB58 on the basis of millimeter and submillimeter observations (Baker et al. 2001 



Sawicki 20011 ). 



Reddv et al. f20ld . 



200a) among UV- 



A similar discrepancy is observed by 
selected galaxies with best-fit stellar population ages of t < 100 Myr. As shown 
in Figure [13] (left), the measured ratio of Lfir/Luv for these young systems falls 
significantly below the local starburst relation, given their observed rest-frame UV 
slopes, p. Both cB58 and the Cosmic Eye are described by similarly young ages as 
well. Therefore, it appears that these young systems are described by a different 
extinction law, one that is steeper (as in the case of the SMC curve), such that a 
given observed amount of reddening corresponds to less attenuation in the rest- 
frame UV. More work is needed to develop a self-consistent evolutionary scenario 
to explain the reddening and attenuation properties of star-forming galaxies as 
a function of the maturity of their stellar populations. Finally, we call attention 
to recent results b ased on Herschel PA CS 160^m measurements of star-forming 



galaxies at z ~ 2 (JNordon et al. 20101 ). A comparison of PACS and UV, dust- 
corrected star-formation rates indicates that the UV-corrected values are, on 
average, 0.3 dex higher (with a scatter of 0.35 dex). More extensive comparisons 
with upcomin g PACS far-IR measurements will be vital for further testing of the 



Calzetti et al. 



law. 
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Related to the overall reliability of the iCalzetti et al.i law, we must also con- 
sider the question of differential extinction of the radiation from stars and ion- 
ized gas. In local starbursts, emission-line tracers of ionized gas appear to 
be systematically more attenuated than the stellar continuum at similar wave- 
lengths. The relationship derived betw een the stellar and nebular extinction is 



E{B - V) stars = 0.4:4E{B - V)nebuiar dCalzetti et al. 200m . Currently, there is 



conflicting evidence about the relative extinction of stars and gas at z ~ 2. Com- 
paring star-formation rates inferred from rest-frame UV an d Ha luminosities 
for a sample of z ~ 2 UV-selected galaxies (see Section 15. 4p , lErb et al. f2006bl) 



conclude that E{B — V)stars ~ E{B — V)nebuiar, and that a .Calzetti et al. 



law applied to correct both UV-continuum and Ha measurements gives rise 
to th e best agreement between sta r-formation rate indicators. On the other 



hand, 



Forster Schreiber et al. (20091 ) find for a set of star-forming galaxies with 



VLT/SINFONI integ ral-field u nit (IFU) maps of Ha emission at roughly the 

I 

same redshift as the 



Erb et al.l sample (see Section 16.2. ip that the best agree- 
ment between Ha and UV-derived star-formation rates results when Ha lumi- 
nosities are corrected by an additional factor of -^ 2, and with the assumption 
E{B — V) stars = 0A4:E{B — V)nebuiar- To Settle the question of differential ex- 
tinction, much larger samples of objects are required with both measurements 
of multiple Balmer emission lines and rest-frame UV estimates of star-formation 
rates. Assembling these measurements will be possible with the next generation 
of multi-object near-IR spectrographs presently coming on-line on 8 — 10-meter 
class telescopes. 
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5.2 Dust Emission 

In addition to absorbing ultraviolet and optical radiation from stars, dust re-emits 
at IR and submillimeter wavelengths. Direct observations of this re-radiated 
emission at long wavelengths have opened a window into the nature of dust 
in distant galaxies, in terms of its temperature, composition, and the sources 
heating it. At far-IR and submillimeter wavelengths, direct observations of dust 
emission are restricted to the most luminous sources, with L^oi > 10^'^ Lq. While 
mid-IR imaging observations have been obtained for lower-luminosity systems 
using Spitzer/MIPS, spectroscopy has been limited to the most luminous sources, 
selected on the very basis of their bright submillime ter or mid-IR emis sion [except 



in cases of strongly gravitationally-lensed systems (jSiana et al. 20091 )]. 

Until recently, direct measurements of dust temperatures in z > 2 ULIRGs 
only existed for small samples. Indirect estimates of the dust temperature, T^, in 
SMGs were obtained by measuring the flux at individual rest-frame far-IR (i.e., 
observed 850 /im) and radio (i.e., 1.4 GHz) wavelengths, and ass uming that th e 



local correlation between far-IR and radio luminosities applies ([Condon 19921 ). 
Based on the ratio between 850 /im flux and inferred total far-IR luminosity, T^i 
can be inferred, assuming that the dust emission follows a single-temperature 

with the emis- 



modified blackbody spectrum of the form Sy oc c:im(hL -,'hT ) i ' 



sivity, 13, set to a value of 1.5. 



Chapman et al. (20051 ) use such a method to 



characterize their sample of 73 SMGs with spectroscopic redshifts at a median 
redshift of z = 2.2. For this sample, the median inferred dust temperature is 
Td = 36-o7K, ~ 5 K cooler than local ULIRGs with similar IR luminosities. In 
order to constrain the dust temperature more directly, photometric measurements 
at multiple rest-frame far-IR wavelengths are required, and the observations be- 
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tween Spitzer/MIPS and SCUBA/850 /im wavelengths have until recently been 

limited. These include Caltech Submillimeter Observatory (CSO) Submillime- 

ter High Angular R esolut ion Camera (SHAR C-2) observations at 350 /xm from 



Kovacs et al. (20061 ) and 



CoDPin et al. (2008i ) for a total of ~ 30 SMGs with 



previous 850/im detections. In these studies the additional far-IR SED point 
suggests a median T^ = 30 — 35 K, consistent with the earlier, indirect estimate 
of Trf = 36 K for SMGs. Furthermore, the total far-IR luminosities of SMGs in 
these samples are better constrained and typically LpiR ~ few x IO^'^^Lq, with 
dust masses MduM ^ lO^M^?^, signific antly larger than those observed in local 



starburst galaxies (jCoppin et al. 20081 ) 



New Herschel observations at 100 and 160/im with PACS and at 250, 350, 
and 500/um with the Spectral and Photometric Imaging Receiver (SPIRE) have 
recently provided significantly more refined estimates of dust temperatures in 
high-redshift U LIRGs, and revealed a temperature diver sity hinted at by earlier. 



indirect studies (jCasev et al. 2009i . iChapman et al. 20041 ) . 



Combining data in the 



Chapman et al. (20ld ) 



three SPIRE channels with 850/im SCUBA measurements 
find a median of T^ = 34 it 5 K for a sample of 31 SMGs and a hotter median 
r^i = 41 lb 5 K for 37 radio-selected ULIRGs with fainter submillimeter fluxes 
(where the errors represent the standard deviations of the samples, not of the 
median values). These samples are characterized by median far-IR luminosities 
of Ljfi = 7.1 X IO^^Lq and 3.8 x W^^Lq, respectively. Such large luminosities 
were inferred previously on the basis of much more limited rest-frame far-IR 
data, with large uncertainties on the total luminosity due to its strong depen- 
dence on assumed dust temperature. The constraints on the far-IR SED shape 
now offer a much more robust indication of the dust luminosity and temperature. 
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Magdis et al. (2010h ) focus on combined PACS and SPIRE observations for a 



sample of star-forming MIPS 24/xm-selected ULIRGs at z ~ 2, finding a broad 
range of dust temperatures spanning from 25 < T^^ < 65 K, with a median of 
Td = 42 K, and a range of IR luminosities of 1 TxlO^^Lfr, < Ljr < 8 .7 x 10^^^^ 



Chapman et al. (20041 ) , it is shown 



0- 



In particular, and as emphasized previously by[| 
that ULIRGs with hotter dust temperatures and fainter submillimeter fluxes at 
850 — 1200/im would be missed in current ground-based surveys of SMGs, given 
their sensitivity limits at ^ 1mm. Therefore, ground-based surveys for z > 2 
ULIRGs selected on the basis of 850//m fluxes may be biased towards the coolest 
of the most luminous galaxies. The origin of the diversity in dust temperatures 
among these most luminous systems is an open question. 

While the far-IR SED offers constraints on the thermal properties and total 
luminosities of the most luminous galaxies at z > 2, the mid-IR spectral range 
reveals the nature of their smaller dust grains and PAH molecules. Furthermore, 
a fundamental issue regarding ULIRGs at high redshift concerns the nature of 
the sources heating the dust that emits such copious amounts of far-IR radiation. 
The two basic alternativ es are star-formation or AGN activity. Whil e rest-frame 



UV and opti cal spectra (IChapman et al. 2005 



Swinbank et al. 20041 ) and X-ray 



observations ([Alexander et al. 20051 ) of SMGs indicate evidence for AGN activity 
in these systems, it is challenging to assess the bolometric importance of the 
AGNs from these data. The rest-frame mid-IR spectral region, from ~ 5 — 15/xm, 
offers strong discriminatory power between star-formation and AGN activity, and 
constraints on the relative contributions of each to the bolometric luminosity of 
high-redshift ULIRGs. The energetics of the radiation field are determined by 
the prominence of PAH emission features at 6.2, 7.7, 8.6, and 11.3/im, relative 
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to the strength of an underlying power-law continuum from hot dust and silicate 

dust absorption at 9.7/im. The PAH emission is primarily tied to star formation, 

whereas the hot dust continuum is associated with AGN activity. The strength 

of silicate a bsorption potent i ally indicates the importance of a buried nuclear 



component (jPope et al. 20081 . 



Saiina et al. 20071 ). 



The IRS instrument onboard Spitzer has proven critical for untangling the 
processes powering the dust emission in ULIRGs at z ~ 2. Strikingly, the 
mid-IR spectra of the majority of SMGs exhibit s trong PAH emission features 



( Menendez-Delmestre et al. 2009 



Pope et al. 20081 ). indicating the dominance of 



star-formation over AGN activity in these systems. While many SMGs have X-ray 
properties suggesting the presence of an AGN, this component does not appear 
to be energetically dominant in terms of the bol ometric luminosity, typically con- 



tributing < 30% of the luminosity in the mid-IR (JMenendez-Delmestre et al. 20091. 



Pope et al. 20081 ) . Only 4 out of the 24 sources present ed in 



and 2 out of 13 of those described in lPope et al. (20081 ). have mid-IR spectra dom 



Menendez-Delmestre et al. (200^ 



inated by a power-law continuum, and the composite spectra for both samples 
exhibit pronounced P AH features. On th e other hand, the z ~ 2 24/xm-selected 



ULIRGs analyzed by 



Saiina et al. (20071 ) are more heterogeneous in the mid- 



IR. The majority (~ 75%) of these sources have mid-IR spectra dominated by 
a power-law continuum and therefore AGN activity. At the same time, more 
than half of these power-law sources have PAH emission features as well, indi- 
cating contributions from both AGN and star-formation processes. The minority 
(~ 25%) of PAH-dominated sources also indicate evidence for AGN activity in 
the form of hot dust continuum. Furthermore, the significant strength of 9.7/im 
silicate absorption in ~ 25% of the sample indicates the presence of an obscured. 
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compact nuclear c omponent. 



Figure [T^ from lPope et al. (20081 ) , illustrates the diversity among high-redshift 
ULIRG mid-IR spectra, including the PAH-dominated flavor common among 
SMGs and in a minority of 24/im-selected ULIRGs, and the power-law domi- 
nated spectra common in 24;um-selected ULIRGs, with different amounts of sili- 
cate absorption. While SMGs and 24/im-selected ULIRGs appear to have similar 
bolometric luminosities, AGNs appear to play a more significant role in the 24/im- 
selected ULIRGs, which also have hotter dust temperatures. These differences 
have been interpreted in terms of an evolutionary scenario, in which SMGs and 
24/im-selected ULIRGs re present, respectively, earlier and 



rich, major merger event (jPope et al. 20081 . lYan et al. 2010l ). Through the pro- 



ater stages of a gas- 



gression of these stages, obscured nuclear AGN activity grows in importance as 
the system evolves into an unobscured QSO, and, eventually, a massive elliptical 
galaxy. While this proposed scenario is intriguing, additional constraints on the 
number densities, stellar mass and gas content of these various high-redshift sarn - 



ples is necessary to establish robust connections between them (JYan et al. 2010l ). 



5.3 Molecular Gas Content 

For a complete characterization of the process of star formation at high redshift, 
observations of the molecular phase of the ISM are critical. Stars form directly 
from this dense interstellar phase, which dominates the cool gas content of the 
most actively star- forming galaxies f|Blitz &: Rosolowskv 20061 ). In the local Uni- 
verse, the star-formation rate per unit area, T,sfr, is directly related to the total 
gas surface density, T,gas , according to the empirical Schmidt L aw, T^sfr ot S^ 



with observational determinations of N ranging from 0.9 — 1.7 (JKennicutt 19981 ) 
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Tracing this relationship at high redshift is key to understanding how gas is con- 
verted into stars during the epoch of peak star formation. Furthermore, the 
star-formation rate and gas content can be related to characterize the efficiency 
of star formation, as well as the timescale on which gas will be depleted. With 
short (compared to the Hubble time) gas depletion timescales, star formation 
can only be sustained by the ongoing accretion of gas from the IGM. The rela- 
tionship between star formation and the available gas reservoir, in terms of its 
mass, baryonic mass fraction, and spatial extent, therefore offers crucial inputs 
into models of galaxy formation, which must include a proper description of the 
balance between gas accretion, the conversion of gas into stars, and the energetic 
feedback related to the process of star formation. 

In the local Universe, rotational transitio ns of carbon monoxide (CO) serve as 



excellent tracers of molecular hydrogen gas ([Young fc Scoville 19911 ) . These same 
transitions are used to trace the molecular gas content of high-redshift galaxies, 
using millimeter- wave and radio telescopes. The first such observations were of 
extreme sources with known ultraluminous far-IR luminosities , such as QSQs and 



SMGs, and were su ggestive of large amounts of dust and gas (JGenzel et al. 2003 



Qmont et al. 19961 ). Recently, however, CO detections have been achieved for less 
extreme systems, well into the LIRG regime and falling on the correlation between 
SFR and Mstar that appears to describe more quiescently st ar-forming galaxies 



over the range lOM rr^yr'^ < SFR <~ few x lOOMeyr-^ (IDaddi et al. 20071 . 



Tacconi et al. 20081 ). For observations of CO at high redshift, the IRAM Plateau 
de Bure Interferometer (PdBI) has played a dominant role, tuned to detect- 
ing various upper-level transitions [e.g. C0(2-l), CO(3-2), CO(4-3), and so on] 
in the millimeter range of the spectrum. Very recently, longer-wavelength ob- 
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servations carried out at the Expanded Very Large Array (EVLA) and Green 

Bank Telescope (GBT) have been tuned to the ground-sta te CQ(l-O) transi- 



tion, which is mo re directly tied to the molecular gas mass f|Harris et al. 201C1 . 



Ivison et al. 2011 



In order to infer the mass of molecular gas, Mgas, associated with detected CO 
emission, there are two major sources of systematic uncertainty. One is in the con- 
version factor between CO(l-O) emission luminosity and Mgas- The second results 
from the fact that, for the most part, high-redshift CO observations have been of 
upper-level J transitions, whereas the CO-to-i/2 conversion factor is calibrated 
for the CO(l-O) transition. The ratios between upper-level and ground-state J 
transitions depend on the excitation and physical conditions in the molecular 
gas, and incorrect assumptions about these conditions will result in a bias in the 
inferred CO(l-O) luminosity. In some cases, higher-J transitions, which are sen- 
sitive to warmer and denser regions, may even offer an inco rrect representatio n 



of the spatial distribution of the full molecular gas reservoir ^Ivison et al. 201 ll ) . 



As for the C0-io-H2 conversion factor, a, it has been calibrated in the Milky 
Way with a value of umw ~ 4 — 5, in units of Mq{K km s~^ pc^)~^. This value 
is very close to the theoretical expectation based on the assumption that CO line 
emission is produced in discrete, vir ialized clouds obeying sc aling relations among 



their masses, sizes, and linewidths (JYoung &: Scoville 1991 



). On the other hand. 



in galactic nuclei and local starburst galaxies, different dynamical and geometri- 
cal conditions apply to the molecular gas, which may reside in a smoother, diskier 
configuration whose motions are additionally affected by the gravitational poten- 
tial from stars. These differences result in a lower conversion factor for starbursts 



of a = 0.8 — 1.6 (JTacconi et al. 20081 ). For observations of high-redshift SMGs, 
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the starburst conversion factor has been adopted (jTacconi et al. 2008. . i2006. ) . In 

fact, using the higher, Milky Way conversion factor would result in baryonic (gas 

plus stellar) masses in excess of the measured dynamical masses. On the other 

hand, for the lower-luminosity LIRGs (both UV-selected and BzK sources), a 

Milky W ay type conversion factor is favored on the basis of similar dynamical ar- 



guments (JDaddi et al. 201Q ) and the idea that CO emission in these LIRGs arises 



in viriali zed clouds with dens ities similar to those observed in local quiescent disk 



galaxies (jTacconi et al. 20081 ) 



We now consider the excitation of the molecular gas, which will affect the 
conversion between higher-J transitions and CO(l-O). While the CO levels in 
SMGs appear to be thermally populated up to J > 3, simultaneous observations 
of CO(l-O) (VLA), and C0(2-l), and CO(3-2) (PdBI) transitions in a BzK 
LIRG at z ~ 1.5 suggests that the CO(3-2) level is significantly subthermally 
excited, similar to what is observed in the Milky Way and other local disk galaxies 



(JDannerbauer et al. 20091 ). These low-excitation physical conditions require a 



larger conversion factor from CO(3-2) to CO(l-O), and therefore, a larger inferred 
molecular gas mass for a given CO (3-2) line luminosity. 

CO observations of both ULIRGs and LIRGs at z ~ 2 have yielded many im- 
portant insights into the nature of star formation in different regimes of total 
galaxy luminosity. Before even considering inferred molecular gas masses, it is 
worth emphasizing the relationship bet ween far-IR luminos ity, LpiR, and CO lu- 
minosity, Leo- As shown in Figure [TSllGenzel et al. (2010l ) have assembled mea- 



surements of local star-forming galaxies (both quiescent and merging systems), 
as well as z ~ 1.5 — 2.0 UV-selected ("BX") and BzK sy stems and (more IR - 



luminous) SMGs. Most strikingly (and also pointed out by 



Daddi et al. (20ld )^ 
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at a fixed Leo, high-redshift SMGs appear to produce 4 — 10 times more far-IR 

luminosity than their UV-selected and BzK counterparts, which cannot simply 

be attributed to enhanced AGN activity in SMGs. This trend mirrors the one 

observed between local mergers and more quiescently star-forming galaxies. 

In terms of star-formation rate and molecular gas surface densities, SMGs are 

characterized by significantly higher T,sfr at a given Y^gas- The UV-selected 

and BzK galaxies follow relations between T,sfr and 'Sgas that are similar to 

the one observed among quiesc ently star-forming g alaxies in the local Universe, 



with a slope of A^ = 1.1 — 1.2 (jGenzel et al. 20ld ). SMGs follow an analogous 



relation but with a hig her overall normaliz ation, perhaps reflective of their shorter 



dynamical timescales (JGenzel et al. 2010 ) , and also of processes that increase the 



efficiency of star forma tion in the turbule nt, merger-driven environment that may 



be common in SMGs (JEngel et al. 20101 ). More speculatively, a top-heavy IMF 



in SMGs may produce more far-IR luminosity for a given mass of stars formed. 
However, the evidence for IMF variat ions at high-redshift is only in direct at 



this point, and not specific to SMGs (JDave 2008 



same time. 



van Dokkum 2008 ). At the 



Ivison et al. (20111 ) offer a potential caveat regarding the inferred 



bimodality of star-formation efficiencies among SMGs and other, more quiescent 
z ~ 2 systems. In so far as the Leo values are based on a range of CO transitions, 
the conversion of these to CO(l-O) relies on a proper characterization of the 
gas excitation conditions. Larger samples of z ~ 2 star-forming galaxies with 
uniform CO(l-O) measurements will be crucial for characterizing the diversity 
(and actual bimodality) among star-formation efficiencies in high-redshift star- 
forming systems. 

The inferred molecular gas reservoirs in UV-selected and BzK galaxies are 
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extended on scale s of several kpc, wit h typical masses of 0.5 — 1.0 x 10^^ Mq 



(Daddi et al. 201C 



Tacconi et al. 20 id ). These correspond to median gas frac- 



tions of 0.44 and 0.57, respectively, for the 10 UV -selected and 6 Bz K galaxies 



and typical gas depletion timescales of ~ 0.5 Gyr (JGenzel et al. 2010l ). The typi- 
cal size of the molecular gas distribution is smaller (half-light radii of ~ 2 — 3 kpc) 
in SMGs, and, while the molecular gas masses and fractions are similar on av- 
erage, the apparent gas depletion scales are a factor of several shorter than for 
the UV-selected and BzK galaxies, due to the higher star-formation rates. With 
the Atacama Large Millimeter Array (ALMA) it will be possible to extend these 
studies down to fainter luminosities and characterize the gas reservoirs of more 
typical systems at z > 2 that make up the bulk of the star formation at those 
early epochs. 

5.4 Metal Content 

While molecular gas provides the material out of which stars form, heavy ele- 
ments constitute an important product of star formation, returned to the ISM by 
supernova explosions and stellar winds. As such, the metal content of galaxies 
reflects the past integral of star formation, modified by the effects of gas inflow 
(i.e. gas accretion) and outflow (i.e. feedback from star formation or black-hole 
accretion). The relative abundances of different chemical elements also provides 
clues about the past history of star formation, in so far as so-called a elements 
(e.g., O, S, Si, Mg) are produced primarily in Type II supernovae events, on 
short timescales (~ 10 Myr), while Fe-peak elements (e.g., Fe, Mn, Ni) are pro- 
duced primarily in Type la supernovae over longer timescales (~ 1 Gyr). The 
metal content of galaxies is especially meaningful when considered in concert with 
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their stellar and gas masses, since the relationships among these quantities ~ and 

deviations from "closed-box" expectations provide constraints on the nature of 

large-scale gas flows (in both directions). 

There are many different methods for measuring the metallicities of galaxies 
at high redshift, probing both their stellar and gaseous components using rest- 
frame UV and optical spectroscopic features. Stellar metallicity is measured from 
absorption lines, while the metal content of the interstellar gas can be gauged 
from either absorption lines arising in the neutral and ionized ISM, or emission 
features originating in H II regions. 

In general, the continuum S/N and spectral resolution obtained for the rest- 
frame UV and optical spectra of typical z > 2 galaxies are not sufficient for 
robust absorption-line metallicity measurements in individual objects (even us- 
ing 8 — 10-meter class telescopes). The rest-frame UV spectra of star- forming 
galaxies include a host of interstellar features arising from neutral hydrogen and 
both neutral and ionized metal species, but only the strongest, highly-saturated 
interstellar absorption lines are detected on an individual-object basis within 
high-redshift samples (jShaplev et al. 20031 ). These saturated features are not 
useful for metallicity estimates. In exceptional cases of the spectra of strongly 
gravitationally-lensed objects, for which both the continuum S/N and resolu- 
tion are at least an order of magnitude better than average, weak, unsaturated 
interstellar metal absorption features are detecte d and can be used for inter- 



stellar metallicity estimates. 



Pettini et al. (20021 ) measure weak features from 



a, Fe-peak, and intermediate elements (i.e. nitrogen) in the spectrum of the 
gravitationally-lensed z = 2.73 galaxy, cB58 (see Section [5. ip . to infer an actual 
interstellar abundance pattern. Based on the relative enhancement of a to both 
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Fe-peak elements and nitrogen, IPettini et al. (2002, ) estimate a "young" age for 
cB58, of less than ~ 300 Myr, the timescale for nitrogen enrichment - an unusual 
case of a chemical constraint on the past history of star formation. 

Rest-frame UV stellar absorption features from both the photospheric and 
wind features o f hot stars can in principle be used to infer stellar metallicity. 



Rix et al. (20041 ) develop calibrations for metal-sensitive stellar photospheric ab- 
sorption indices at 1370, 1425, and 1978 A. These have been applied to a small 
number of individual lensed an d unlensed star-forming galaxies at z '^ 2 — 3 



(jDessauges-Zavadskv et al. 2010l . 



Quider et al. 2009 



Steidel et al. 2004l l. yield- 



ing results consistent for the most part with other metallicity indicators. On 
the other hand, the 1978 A Fe III index was measured in a composite spectrum 



of 75 star- forming galaxies drawn from the GMASS survey (JHallidav et al. 20081 ). 
in fact suggesting a systematic enhancement of a relative to Fe, when compared 
with the expected oxygen abundance for objects of the same stellar mass (see 
below). The shape of the C IVA1549 P-Cygni wind feature from O and B stars 
is also sensitive to metallicity (as well as the form of the IMF), and has also 
been used to estimate stellar metallicity in a few (mainly gra yitationally- lensed) 
^; ^ 2 — 3 galaxies with a dequate S/N and spectral resolution (jPettini et al. 200d . 



Quider et al. 2009 



20ld ). These interstellar and stellar absorption metallicities 
offer intriguing and detailed probes of small numbers of special high-redshift 
galaxies, but await the power of future 30-meter-class ground-based telescopes 
for application to large samples of individual, unlensed objects. 

Most results about the metal content of high-redshift galaxies are based on 
measurements of rest-frame optical emission lines from H II regions. These in- 
clude combinations of hydrogen recombination lines (Ha, II/3), and collisionally 
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excited forbidden lines from heavy elements such as oxygen ([OIII], [Oil]), ni- 
trogen ([Nil]), and neon ([Nelll]). At z > 2, rest-frame optical features shift 
out of the observed optical range and require near-IR spectroscopic observations. 
The study of rest-frame optical emission from H II regions at high redshift to 
date has relied on both long-slit and IFU spectroscopy, using instruments such 
as NIRSPEC and OSIRIS at the Keck Observatory, ISAAC and SINFONI at 
the VLT, MOIRCS on Subaru, and GNIRS on Gemini-South. Emission lines 
are primarily used to infer the gas-phase abundance of oxygen (expressed a s 



12 + log(0/H), where the solar value in these units is 8.66; 



Asplund et al. 2004 ) 



and are based on the relations between particular sets of emission-line ratios 



and metallicity. T hese relations have been both emp i rically calibrated in t 



re lo- 



cal Universe (e.g., 



Nagao. Maiolino Sz Marconi 2006 



Pettini k Paeel 2004'! and 



theoretically modelle d using photoionization codes (e.g.. 



Kewlev k Dopita 20021 . 



Tremonti et al. 2004 ). Two commonly-used indicators of 12 + log(0/H) at high 
redshift are the so-called A^2 index, defined as log([NII]A6584/IIa), and R23 = 
([OIII] + [0II])/H/3. From the ground, R23 can be measured within near-IR win- 
dows of atmospheric transmision for various redshift intervals between z ~ 2 and 
z ~ 4, while A^2 is only measurable up to z ~ 2.6, at which point the thermal 
background becomes preventively high. 

In addition to the significant scatter among different emissio n-line indicators 



(differences as large as 0.7 dex in metallicity for the same galaxy) (JKewlev fc Ellison 20081 ) 
another potential source of bias when using locally-calibrated metallicity indica- 
tors to interpret the emission-line ratios of high-redshift galaxies is the assumption 
that the physical conditions in high-redshift galaxy H II regions are similar to 
those in local galaxies. Small samples of UV-selected galaxies at z ~ 2 with 
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measurements of both [0III]/H/3 and [NII]/Ha, and/or [OIII]/[OII] and R23, 



i ndicate systemat i c offsets from the ex citation sequence of low-redshift galaxies 

J) 



(Erb et al. 2006a 



Hainhne et al. 20091 ). These differences may be indicative of 



a systematically higher ionization parameter and/or electron density, or harder 
ionizing radiation spectrum, with correspondingly different translations between 
empirical line ratios and physical metallicities. Larger samples of objects at z ~ 2 
with such measurements will be required to determine the origin of these apparent 
offsets, relative to local star-forming galaxies. 

At this point, N2 measurements have been obtained for a wide variety of z ~ 2 



sourc es, numbering 



Ce-g-. 



Erbet al. 2006a . 



Yoshikawa et al. 20ld ) to SMGs ( 



200 and ranging from UV- and near- 



Havashi et al. 2003. 



Kriek et al. 2007 



Swinbank et al. 2004 ). The samp 



R selected galaxies 



Shaplev et al. 2004 



with emission-line metallicity measureme nts at z > 3 is much smaller ([Maiolino et al. 20081 . 



Mannucci et al. 2009 



Pettini et al. 2001 



e of objects 



) , with metallicities based on R23 or other 



combinations of oxygen, [Nelll], and Balmer lines for ~ 30 objects. 

The measurement of galaxy metallicities along with stellar masses allows for 
the construction of the galaxy mass-metallicity (or Mgtar-'Z') relation. This re- 
lations hip has been studied in the local Universe for > 50, 000 galaxies with 
SDSS (jTremonti et al. 20041 ). and used as a tool to constrain the importance 
of inflows and outflows as a function of galaxy mass. Another necessary com- 
ponent for this type of analysis is an estimate of the gas mass, which enables 
a calculation of the metallicity as a function of gas fraction, or, equivalently, 
the effective yield, ye//- Given the small sample of galaxies with direct atomic 
and/or molecular gas measurements both locally and at high redshift, gas masses 
have typically been estimated indirectly on the basis of T,sfr values and galaxy 
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sizes, and an assumption that the locally-calibrated relation between Y^sfr and 



Sgas applies (JErb et al. 2006a 



Tremonti et al. 2004) . We note here that the av- 



erage gas fraction of ~ 50% inferred by lErb et al. (2006d ) for UV-selected star- 



forming galaxies at z ~ 2 agrees well with the CO-based estimate of 44% from 



Tacconi et al. f20ld ). 



The largest surve y of rest-frame optical emission lines at z > 2 is presented 



m 



Erb et al. (2006al ). where spectra covering the Ha region were obtained for 87 
UV-selected galaxies. As shown in Figure [T6t the empirically-measured N2 indi- 
cator increases monotonically in bins of increasing st ellar mass, corresp onding to 



an increase in 12-|-log(0/H). At a given stellar mass. iErb et al. (2006al ) find that 



z ^ 2 galaxies are 0.3 dex lower in metallicity than the local galaxies studied in 



Tremonti et al. (20041 ). Based on indirect estimates of gas masses and gas frac- 
tions, these authors compare their observations with simple chemical evolution 
models. The shallow slope of the z ~ 2 Mgtar — Z and gas-fraction— Z relations 
are then used to demonstrate the importance of galax y-scale wi nds with mass- 



Erb_200a). The observed 



outflow rates roughly equal to the star-formation rate (y 

Mstar — Z slope at z ~ 2 has been viewed as e vidence in favor of a "rn omentum- 

driven" feedback recipe in the simulations of iFinlator &: Dave (20081 ). however 



we caution against placing too much significance on the precise observed value of 
this slope, given the known limitations of the A'^2 indicator (which saturates at 
ro ughly solar metallicity ). 



Maiolino et al. (20081 ) construct the Mgtar — Z relation for 9 galaxies at z 



3.5, showing that, at fixed stellar mass, galaxies at z ^ 3.5 are 0.4 dex lower 



in 12 + log(0/H) than at z ~ 2. 



Mannucci et al. (2010l ) explain the observed 



evolution of the Mgtar — Z relation up to z ~ 2.5 in terms of galaxies populating 
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a more general Mstar — Z—SFR relation, with no redshift evolution. Ac cording to 

this "Fundamental Metallicity Relation" (FMR; iMannucci et al. 20101 ) . galaxies 



at fixed stellar mass but higher star-formation rate will have lower metallicity, 
resulting from the interplay between infalling IGM gas and outflowing enriched 
material. The rapid apparent evolution in the Mgtar — Z relation beyond z ^ 2.5 
requires some evolution in the FMR, yet the sample of galaxies at these redshifts 
with metallicity measurements is too small to draw definitive conclusions. 

In order to use the Mgtar — Z relation to constrain models of gas inflow and 
outflow in star-forming galaxies at z ~ 2 — 4, signiflcantly larger (i.e. at least 
an order of magnitude) samples of galaxies with metallicity measurements are 
required, which are complete down to a given stellar mass or star-formation 
rate. Furthermore, constructing the relations between metal and baryonic (stel- 
lar and gas) content for galaxies with direct estimates of molecular gas con- 
tent will provide more robust estimates of the metallicity as a function of gas 
fraction. Spatially-resolved estimates of chemical abundance gradients will offer 
additional constraints on the inflow and outflow of gas, and the build-up of the 
galaxy stellar population. These have been now reported for a handful of objects , 



with conflicting r esults about the sign of the radial gradient (jCresci et al. 2010l . 



Jones et al. 20101 ). Deeper data for larger samples will be required to settle this 



discrepancy. The next generation of near-IR multi-slit spectrographs and IFUs 
should easily enable these observations. 

6 STRUCTURAL AND DYNAMICAL PROPERTIES 

The diversity among galaxy structural and dynamical parameters in the local 
Universe reflects a range of mass assembly histories, and the relative importance 
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of mergers and smoother mass accretion. Furthermore, properties such as size, 

concentration, mass-surface density, and velocity dispersion are correlated with 

indicators of the nature of stella r populations, such as luminosity, stellar mass, 



and specific star-formation rate (JKauffmann et al. 2003bl . 



Shen et al. 2003l l. Re- 



producing the observed connections among these galaxy structural and stellar 
properties and their redshift evolution within a unified framework constitutes 
a crucial test for models of galaxy formation. In this section we review recent 
results about the structural and dynamical properties of high-redshift galaxies. 
Due to the small apparent sizes of most z >2 galaxies (< 1"), i/S'T imaging and 
ground-based IFU maps assisted by adaptive optics (AO) are both critical for 
resolving overall shapes and ~kpc-scale fine structure. Interferometric observa- 
tions of CO line emission have also provided valuable insight into the dynamical 
properties of the most extreme sources at high-redshift, such as SMGs. While 
high-resolution imaging and dynamical maps are e ach powerful probes individu - 



ally, combining these two methods of observations (jForster Schreiber et al. 2011 



will provide the deepest insights into the nature of the galaxy assembly process 
at high redshift. 

6.1 Structural Properties 

In general, the morphologies of galaxies at z > 2 are characterized by two im- 
portant differences with respect to those of local galaxies. First, the traditional 
Hubble sequence of regular spirals and elliptical galaxies has not settled into place 
by z ~ 2, and a much higher freque ncy of clumpy, irregular morphologies is ob - 



served among star-forming systems (JLotz et al. 2006 



Ravindranath et al. 20061 ) 



Second, both star-forming and quiescent galaxies are more compact at fixed stel- 
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lar mass or rest-frame optical luminosity (jBuitrago et al. 2008 ) . In order to 
characterize the structural properties of high-redshift galaxies, both paramet- 
ric and non-parametric methods have been employed. In the former, model 
Sersic profiles are fit to the data, yielding best-fit values for n, the Sersic power- 
law index, and r^, the effective radius, within which half the galaxy luminos- 
ity is emitted. Sersic profile fits work best for regular, axially-symmetric mor- 
phologies with well-defined centers, so, in order to address the clumpy, irreg- 
ular structures observed among high-redshift systems non-parametric statistics 
such as the Gini, M2n, Multiplicity, and CAS coefficients have been developec 



(lAbraham. van den Bergh &: Nair 20031: 



Conselice. Chapman &: Windhorst 20031 : 



Law et al. 2007 



Lotz et al. 20061 : 



Lotz. Primack &: Madau 2004 ) to characterize 



their concentrations, asymmetries, and dumpiness (i.e. deviations from maxi- 
mally compact configurations). 

The first high- redshift galaxies to be observed with HST were UV-selecte d 



galaxies at z ~ 3 (JGiavalisco. Steidel fc Macchetto 19961 : 



Lowenthal et al. 1997i ). 



The irregular morphologies of these systems, often composed of multiple com- 
pact components and/or irregular nebulosity, were first attributed to the fact 
that the Wide Field Planetary Camera 2 (WFPC2) and ACS are optical im- 
agers, probing rest-frame UV wavelengths at z ~ 3 - a so-called "morphological 
/c-correction." In principle, measurements at these wavelengths might be sensi- 
tive to only the most active regions of star formation rather than the bulk of the 
stellar mass, and potentially affected by patchy dust extinction. However, in prac- 
tice, rest-frame optical imaging with the Near Infrared Camera a nd Multi-Object 



Spectrometer (NICMOS), using both NIC2 and NIC3 cameras dPickinson 200d 



Forster Schreiber et al. 2011 



Kriek et al. 2009 



Papovich et al. 20051 ) and the Wide 
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Field Camera 3 (WFC3) (jQverzier et al. 2010. ) has demonstrated that actively 
star-forming galaxies at z ~ 2 — 3 have very similar rest-frame UV and optical 
morphologies. On the other hand, among the most massive, evolved systems at 
z ~ 2 — 3, there is evidence for more centrally-concentrated, regular morphologies 
at rest-frame o ptical wavelengths, compared with th e structure observed in the 



rest-frame UV (|Cameron et al. 20ld . 



Toft et al. 20051 ). With the recent installa- 



tion of WFC3, it will now be possible to measure rest-frame optical morphologies 
for statistical samples of z ~ 2 — 3 objects at ~kpc-scale resolution. The perfor- 
mance of the IR channel of WFC3 is superior to the NICMOS NIC3 camera in 
terms of resolution, sensitivity, and area. 

While the traditional Hubble sequence of disk and elliptical is not in place 
by z ~ 2, HST imaging has revealed a large diversity of structural parameters 
among galaxies at z ~ 2 — 3, which correlate with their stellar populations in anal- 



ogy with the tre nds observed in t 



he loc al Universe. Bas ed on NICM0S/NIC3 



F160W imaging, IZirm et al. (20071 ) and 



Toft et al. (20071 ) show that quiescent. 



high-redshift z ~ 2 galaxies have systematically smaller sizes and higher stellar- 
mass surface densities t h an actively star- formi ng systems of the sa, me stellar 



mass. 



Franx et al. (20081 ). 



Toft et al. (20091 ). and 



Wilhams et al. (20ld ) use deep. 



ground-based i^-band images to demonstrate the same trend. Based on higher- 
resolution NICM0S/NIC2 F160W imaging for a sample of spectr oscopically- 



Kriek et al. (20091 ) presents a 



confirmed X-selected massive galaxies at z ~ 2, 
comparison of the Sersic profile fits for the quiescent and emission- line objects. 
As shown in Figure [T71 galaxies with redder rest-frame U — B colors, lower spe- 
cific star-formation rates, and SEDs indicative of more evolved stellar popula- 
tions, are characterized by smaller r^ and more concentrated surface-brightness 
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profiles (larger n). This figure also shows the clumpy rest-frame optical mor- 
phologies of the massive, blue galaxies. Objects whose emission-line ratios ac- 
tually suggest AGN activity show stru ctural properties and SEDs th at are sim- 



ilar to those of the quiescent galaxies. 



Forster Schreiber et al. (201ll ) find even 



larger sizes and shallower profiles in NIC2 images of the six clumpy UV-selected 
galaxies in their sample, which also have higher specific st ar-formation rates 

Also similar to the trends observed in the local Universe 

1 

and iForster Schreiber et al. (2011i ) find a correlation between stellar mass sur- 



Franx et al. (;2008l ') 



face density and specific star-formation rate, such that galaxies with higher stella r 



mass densities have lower specific star-formation rates ([Kauffmann et al. 2003bl ). 
We now highlight two important observed morphological phenomena, which 
have sparked much interest from the theoretical community. One is the nature of 
the clumpy morphology and "clumps" in actively star-forming galaxies at z > 2. 
The other is the incredibly compact nature of massive, quiescent galaxies at high 
redshift, and their connection to today's massive, early-type galaxies. 

6.1.1 CLUMPY MORPHOLOGIES IN STAR-FORMING SYSTEMS 
While it is possible to model the rest-frame UV and optical morphologies of z > 2 
galaxies using smooth Sersic profiles, significant residuals result due to the pres- 
ence of non-axisymmetric surface-brightness fluctuations. These "clumps" have 
been fairly ubiquitousl y observed among the rest- frame UV and optical morpholo- 



gies of distant ga . 



Law et al. 2007 



axies (iCowie. Hu fc Songaila 19951 : 



Elmegreen fc Elmegreen 20051 : 



Lotz. Primack &: Madau 20041 ). and are visible among the ac- 



tively star-forming galaxies in Figure [T71 Clumps are characterized by typical 
sizes of ~ 1 kpc, and make up as much as ~ 40% of the rest-frame UV light 



distribution ()Elmegreen et al. 20091 ). While clumpy, irregular structure was first 
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explained as the natural outc ome of the increased prevalence of majo r merger 



events at high redshift (e.g., 



Conselice. Chapman &: Windhorst 20031 ) . various 



lines of evidence suggest t hat other factors may contribute to this phenomenon. 



First, iLaw et al. (20071 ) demonstrated a decoupling between the Gini coef- 
ficient (and any other rest-frame UV non-parametric morphological statistic) 
and star-formation rate. If clumpy morphology was connected with the inci- 
dence of a major merger, the Gini coefficient, sensitive to non-uniformities in 
the light distribution, should displ ay some correlation w ith the parameters de- 



scribing galaxy stellar populations. 



Swinbank et al. (20ld ) also highlight the fact 



that the distribution of rest-frame UV and optical morphologies of SMGs (plan 



sibly the sites of major merger events; 



Engel et al. 201C 



Swinbank et al. 2009 ) 



were statistically indistinguishable from those of UV-selected star-forming galax- 
ies, which are significantly more quiescent in terms of bolometric luminosity and 
star-formation rate. Furthermore, the clumpy phenomenon is far too common 
among star-forming galaxies at z ~ 2 — 3 to be accounted for quantitatively 
by t he predicted rate of major merger eve nts in, e.g., the Millennium Simula- 



tion fjConrov et al. 2008 



Genel et al. 2008 ). Finally, as described below in Sec- 



tion 16.2.11 clumpy structures are commonly observed in galaxies with ordered 
velocity fields indicative of rotating disks - demonstrating that a major merger is 
not necessarily the cause for clumps. Minor-merger interactions and instabilities 
i n gas-rich, turbulent disks are more likely alternatives for the formation of clump s 



(jBournaud fc Elmegreen 2003 : 



Dekel. Sari fc Ceverino 20091 : 



Genzel et al. 20081 '). 



6.1.2 THE COMPACTNESS OF QUIESCENT GALAXIES At fixed 
stellar mass, both star- forming and quiescent galaxies have smaller radii at higher 
redshift. The evolution in the Mgtar — fe relationship is especially dramatic 
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for quiescent systems. As noted by 



Cimatti et al. 2008 



Daddi et al. 20051 . 



many authors (e.g., 



Shaplff 

Buitrago et al. 2008 



Damianov et al. 2009 



Toft et al. 20071 . 



Truiillo et al. 



van Dokkum et al. 20081 ) using rest-frame UV and optical HST imaging and ground- 



based data, z ~ 1.5 — 2 quiescent galaxies have radii that are ~ 2 — 5 times 
smaller than those of local early-type galaxies of the same stellar mass. For ex- 



ample, z ~ 2 quiescent galaxies wit 



h M. 



Tp ~ 



star 



10 M(Ti are measured to have 



1 kpc (jvan Dokkum et al. 20081) . compared wit h ~ 3 kpc in local early- 



type galaxies of the same mass fiDamjanov et al. 2003 ). The observed differences 
in size correspond to factors of ~ 10 — 100 in physical density! Furthermore, 



Taylor et al. (20101 ) demonstrate the extreme scarcity of such compact, early- 
type systems in the local Universe. The resulting challenge from these discover- 
ies consists of relating high-redshift massive, compact quiescent galaxies to their 
low-redshift counterparts - since structural evolution is required for the distant 
objects to resemble the current ones. Scenarios considered to explain the ob- 
served e volution in size include the effects of major and minor diss ipationless 



mergers (jHopkins et al. 201C 



200S 



Naab. Johansson &: Ostriker 20091 ) as well as 



energetic feedback from an accreting black hole (jFan et al. 20081 ). 



Many important caveats have been raised about the observations of apparent 
compactness among quiescent galaxies. These include (a) the possibility of not 
accounting for extended low-surface brightness emission which would tend to in- 
crease the inferred radius; (6) using in some cases rest-frame UV rather than rest- 
frame optical imaging; (c) radially-dependent stellar M/L ratios that increase 
outwards, causing an underestimate of the mass at larger radii; (d) a reliance 



on photometric rather than spectroscopic redshifts; anc 



lar population models used to estimate stellar masses (jvan Dokkum et al. 2008 



(e) uncertainties in stel- 



GaJa xies at z = 2 — 4 
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20ld ). Most of these concerns have been addressed, using samples with spectro- 



scopic redshifts and deeper surface-brightness hmits. Additionahy, crude dynam- 
ical information has been obtained from stacked and individual spectra of com- 
pact, quiescent galaxies, based on an estimate of velocity dispersion from stellar 



absorption lines (jCappellari et al. 2003 : 



Cenarro k Truiillo 2009 : 



van Dokkum. Kriek &: Franx 2009 ) 



The estimated velocity dispersions appear to confirm the large masses inferred 
from stellar population modeling. Furthermore, results at z ~ 1 based on a com- 
parison of radii and dynamical mass estimates for 50 galaxies suggest a factor 
of 2 evolution in size at fixed dynamical mass, c onsistent with the evolu tionary 



trend observed to z ~ 2 based on stellar m asses (ivan der Wei et al. 20081 ') . 



van Dokkum et al. (2010l ) stack the 



Employing a complementary technique, 
ground-based near-IR surface-brightness profiles for rare, massive galaxies of 
constant number density in five redshift bins between z = 2.0 and z = 0.1, 
such that the higher-redshift objects can plausibly represent the progenitors of 
the objects at lower redshift. The observed evolution in these stacked surface 
brightness profiles shows that the mass within an inner region of r = 5 kpc re- 
mains roughly constant, while the mass at larger radii builds up smoothly as 
a function of decreasing redshift. This result suggests the importance of minor 



merger events in the growth of ni assive, elliptical galaxies (JHopkins et al. 200S 



Naab. Johansson &: Ostriker 20091 ). and that the structure of elliptical galaxies 
is not self-similar as a function of redshift. Based on deep near-IR spectroscopy 
with upcoming instruments on the ground and in space, statistical samples of 
dynamical mass estimates for passive galaxies at z > 2 will be crucial for even 
more robustly distinguishing among evolutionary scenarios. 
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Dynamical studies of high-redshift galaxies have mainly been limited to star- 
forming objects, using rest-frame optical emission lines from ionized gas or millimeter- 
wave emission from CO tracing molecular gas. While the sample of existing 
rest-frame UV spectra for h igh-redshift galaxi e s is m uch larger than either of 



these other types of dataset (jSteidel et al. 2003 



2004 ) ■ the emission and absorp- 



tion features typically detected in the rest-frame UV range are broadened by 
non-virial motions from outflowing gas and radiative trans fer effects, and un 
fortunately do not permit useful dyna mical measurements (jPettini et al. 2001 



Shaplev et al. 20031 . 



Steidel et al. 2010l ). The first dynamical probes of high- 



redshift star-forming g alaxies consisted of long-slit near-IR spectroscopy of ~ 100 



UV-selected galaxies (JErb et al. 2003 



2004, 



maps of very small samples of SMGs ( 



2006c. 



Pettini et al. 2001) and CO 



Genzel et al. 2003 



Neri et al. 20031 '). In 



deed, CO observations can be used not only for estimating molecular gas masses 
but also for probing the dynamical properties of a galaxy. These early ob- 
servations revealed a typical Ha or [OIII] emission-line velocity dispersion of 
o" ~ 100 km s~ for UV-selected galaxies. Measured linewidths were ~ 2 — 3 
times larger in the SMG CO maps, which also showed evidence for complex, 
multi-component emission-li ne morphology. A particularly striking result from 



the long-slit Ha studies of 



Erb et al. (2006d ) consists of the measurement of 



spatially-resolved, tilted emission-lines, indicative of velocity shear, and, per- 
haps, rotation - i.e., preliminary evidence for disks at high redshift. Below we 
summarize some of the latest results from both near-IR IFU and CO dynamical 
studies of high-redshift galaxies, highlighting many open questions in this rapidly 
evolving field. 
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6.2.1 INTEGRAL-FIELD UNIT OBSERVATIONS OF TURBULENT 

STAR-FORMING GALAXIES With the advent of both seeing-hmited and 
AO-assisted near-IR IFU spectrographs on Keck, the VLT and Gemini telescopes 
several years ago, the study of high-redshift galaxy dynamics has advanced consid- 
erably. Some striking trends have emerged from a total sample of ~ 100 objects 
with IFU kinematic observations, which are no longer subject to the spatial- 
sampling limitations of long-slit spectroscopy taken at a fixed position angle . We 
focus here on the VLT/SI NFQNI survey of IForster Schreiber et al. (2003 ) and 



the Keck/OSIRIS survey of 



Law et al. (20091 ) , which represent two of the largest 



campaigns at z > 2 using this new observationa l capability. The "SINS Ha Sur 



vey" presented in 



Forster Schreiber et al. (2009 ) includes 62 star-forming galaxies 



at 1.3 < z < 2.6 selected in the rest-frame UV, optical, near-IR, and submillime- 
ter. The sample is representative of massive, Mstar > 10^*^ Mq, star-forming 
galaxies in this redshift range, with a median stellar mass and star-formation 
rate of M,^tnr = 3 x IO^^Mq and SFR = 70 M^yi-^ , respectively. The sample of 



Law et al. (20091 ) includes 12 UV-selected star-forming galaxies at 2.0 < z < 2.5 



with a median stellar mass and star-formation rate of Mstar = 10"*^*^ Mq and 
SFR = 40 M0yr~^, respe ctively. In additio n to having lower stellar masses on 
average, the objects in the iLaw et al. (2009i ) sample are also characterized by a 
smaller median Ha half-light radius than the S INS Ha sample (1.3 versus 3.1 kpc). 



All OSIRIS observations in iLaw et al. (20091 ) are AO-assisted, with an effective 
PSF of ~ 0.15", corresponding to ^ 1.2 kpc. The majority (^ 85%) of the 



SINFONI observations discussed by 



Forster Schreiber et al. (20091 ) are seeing- 



limited, with a typical resolution of ~ 0.6", corresponding to ~ 4.9 kpc; the 
remainder were observed using either laser- or natural-guide-star AO, with reso- 
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lutions ranging from 0.17" — 0.41", corresponding to ~ 1.6 — 3.4 kpc. 

These IFU surveys have revealed a wealth of diversity among the dynami- 
cal properties of distant star-forming galaxies. The key observables are maps 
of velocity, velocity dispersion, and line-intensity, which can then be used as 
inputs to detailed dynamical models, or analyzed in a more empirical sense. 
Figure [18] demonstrates the range of velocity fields observed in the SINS Ha 
survey. Roughly 1/3 of these velocity maps show clear gradients and are clas- 
sified as "rotation-dominated," with corresponding maps of velocity dispersion 
showing central maxima; roughly another 1/3 qualify as "mergers" as evidenced 
by kinemetric modeling (,Shapiro et al. 2009. ) or distinct multiple components; 
finally, 1/3 of the sample are "dispersion-dominated", with no evidence for ro- 
tati on or shear, but characterized by a significant velocity dispersion. Half of 



the iLaw et al. (20091 ) sample belong to this last category, with no evidence for 
rotation. 

Another basic result consists of the high degree of turbulence observed in the 
ISM of high-redshift galaxies, along with a small ratio of rotational to random 
motions. The spatially-resolved velocity dispersion maps indicate local velocity 
dispersion values ranging from aiocai = 30 — 90 km s~^ in t he SINS Ha sample 
and aiocai = 60 — 100 km s^^ in the sample of iLaw et al.l . Even in rotation- 



dominated systems, the velocity dispersion is significant. The median value of 
Vyot/o' is 4.5, which i s considerably lower than the values of ~ 10 — 20 observed in 



local spiral galaxies ([Forster Schreiber et al. 20091 ). Of course, in the dispersion- 



dominated systems, Vrot/<^ is even smaller (i.e., < 1). The prevalence of rotation 
and rotation-dominated systems appears to be correlated with stellar mass, with 
a higher fraction of rotating systems at higher stellar masses. At the same time, 



Galaxies at z = 2 — 4 77 

galaxies \ yith small stellar masses and large gas fraction s tend to have negligible 



rotation (JForster Schreiber et al. 20091 . 



Law et al. 20091 ). 



The star-forming galaxies with IFU observations also appear to follow scaling 



rela tions between their basic dynam ical and stellar properties. Both 



Bouche et al. (200i ) 



and iForster Schreiber et al. (20091 1 present a strong correlation between galaxy 
circular velocity, Vrot, and size for galaxies selected in the rest-frame UV, optical, 
and near-IR. This correlation is indistinguishable from the one describing local 
disk galaxies (despite the significant differences in Vrot/o'), and both UV-selected 
and BzK galaxies follow it. On the other hand, SMGs occupy a distinct re- 
gion of Wr-oi— size parameter space, with significantly larger velocities and smaller 
sizes. This difference potentially reflects the lower angular momenta and higher 
matter densities present in S MGs relative to mo re quiescently star-forming UV- 



selected and BzK systems (jBouche et al. 20071 ). Limit ing their analysis to 18 



Cresci et al. (20091 ) dis- 



rotation-dominated systems from the SINS Ha sample, 
cover a "Tully-Fisher" correlation between stellar mass and Vrot, which has the 
same slope as the correlation among local disk galaxies, but offset by ~ 0.4 dex 
towards lower stellar mass at fixed Vrot- Furthermore, a comparison of dynamical 
and stellar masses indicates a strong correlation, but suggests that galaxies must 
contain significant gas fractions to explain the diff'erences between the two mass 
estimates. In particular, galaxies with young stellar ages appe ar to contain large r 



Erb et al. (2006d ). 



gas fractions, consistent with the previous, long-slit results of 

The spatially-resolved Hq maps indicate that clumpy morphologies are a fea - 



ture of not only continuum surface brightness (lElmegreen &: Elrnegreen 20051) 



but also the emission from ionized gas (jForster Schreiber et al. 2009i . iGenzel et al. 20081 ) . 



Given that star-forming clumps are present even in galaxies with ordered rota- 
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tion fields and no evidence of major merging, and that tlie clumps appear to 

follow the general velocity field of the larger systems in which they are em- 
bedded, it is unlikely that clumps constitute distinct, accreted systems. An 
alternative explanation for observed clump properties is based on predicting 
the size and mass-scale on which gas should fragment in a turbulent, high- 
surface-density disk. Specifically, the high turbulent velocities observed (as de- 
scribed above) lead to an expected Jeans length for fragmentation of ~ 2.5 kpc 



10^ Mq pc" 



(JGenzel et al. 20081'). Folding in t he typical observed disk surface densities (^ 



Genzel et al. (20081 ) estimate clump masses of '^ 10^ Mq, simi- 



lar to what is observed. Numerical sir nulations of turbulent gas-rich disks also 



appear to produce ~kpc-scale clumps (JBournaud &: Elmegreen 2003 ). Theoret- 



ical arguments (both analytical and numerical) suggest that these clumps may 
migrate inwards and coalesce to form a bulge on ^ 0.5 — 1.0 Gyr ti mescales 



(|Dekel. Sari fc Ceverino 20091 : 



Elmegreen. Bournaud &: Elmegreen 20081 ). On the 



other hand, feedback from star formation may disrupt the clumps before they 
reach the center, and they may simply contribute to the overall growth of the 
disk. The role of clumpy structures in the growth of massive galaxies is an ex- 
tremely active area of research; new observational results will help to elucidate 
the origin and fate of the clumps. 

Local velocity dispersions of ~ 50 — 100 km s~ appear to be a generic prop- 
erty of high-redshift star- forming galaxies at z ~ 2. Determining the cause of 
this large apparent turbulence is an important goal for models of galaxy forma- 
tion. Proposed scenarios include the conversion of the gravitational potential 
energy from infalling accreting matter into random kine tic energy or collis ions 



between large clumps as they migrate inwards (see, e.g.. 



Genzel et al. 2011 



, for 



Galaxies at z = 2 — 4 

a review of these possibilities). On the other hand, 



1 
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Lehnert et al. (2009. ) ar- 



gue that mechanical energy feedback associated with star formation - in the 
form of stellar winds and supernovae explosions - is the driving force behind 
the large observed interstellar turbulence. Various authors have compared the 
surface-density of star-formation Tisfr with the local velocity dispersion (pre- 
sumably an indication of the connection between star-formation feedback and 
interstellar turbulence) and have arrived at diffe rent conclusions about what 



the weak observed correlation (if any ) signifies (jForster Schreiber et al. 20091 . 



Genzel et al. 2011 



Lehnert et al. 20091 ). The cause of large turbulent velocities 



in high-redshift star-forming galaxies clearly remains an open question. 
6.2.2 THE NATURE OF SMGS Dynamical information can also be po- 
tentially used to investigate the origin of the extreme luminosities among SMGs. 
Both near-IR IFU maps and CO observations of SMGs have been used in this 
endeavor. The observed CO profiles are broad (FWHM typically several hundred 



km s ^) and often double-peaked ( Greve et al. 20051 ). A broad, double-peaked 



profile can be i ndicative of either a m assive, rotating disk or, alternatively, a 
merger event. 



Swinbank et al. (20061 ) present near-IR IFU observations of 8 
SMGs at 1.3 < z < 2.6, probing rest-frame optical line emission. At least five of 
these systems show evidence for two or more distinct dynamical components, sug- 
gestive of mer ging. Based on su barcsecond-resolution CO maps for 12 SMGs at 
1.2 < z < 3.4, lEngel et al. (201Q ) find evidence in 5 systems for two distinct spa- 
tial components, with mass ratios (when possible to determine) closer than 1 : 3 
— the standard threshold for being considered a major merge r. In the remaining 



systems, morphologies are either disturbed or compact, and 



Engel et al. f201(]| ) 



argue that these represent later-stage, coalesced merger events. A significant frac- 
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tion of SMGs clearly show dynamical evidence for major mergin g. Whether an al - 



Dave et al. 201 



ternative scenario of smooth gas infall and minor mergers (e.g., 
can explain the luminosities of the remaining sources will require more robust ob- 
servations of the space densities, and stellar a nd dynamical masse s of SMGs. The 



extreme matter densities inferred for SMGs (IBouche et al. 20071). their high ap- 



parent star- formation efficiencies (|Daddi et al. 20101 . iGenzel et al. 2010l ) , and the 



deviation of these syst ems from global s caling relations between star-formation 



rate and stellar mass (JDaddi et al. 2007i ). all suggest that they constitute very 



unusual events. It remains an open challenge to relate SMGs to other ULIRGs 
at the same epochs (e.g., the 24^m-selected objects of JYan et al. 20071 ). as well 



plausible descendants among the lower-redshift massive-galaxy population. 

7 CONCLUDING REMARKS 

Our path has traveled through many different wavelength ranges and observa- 
tional techniques along the way to characterizing the stars, dust, gas, heavy 
elements, structure and dynamics of high-redshift galaxies at 2 < z < 4. By 
design, this review has focused primarily on the translation of observed quanti- 
ties to physical ones, instead of making systematic comparisons with particular 
theories (although at times, when appropriate, a connection was made). At the 
same time, the rapid development of observations and their physical interpreta- 
tion over the previous decade has yielded an important body of data for input 
into models of galaxy formation (i.e., the diversity of star- formation histories 
and structures among the highest-mass systems; the connection between stel- 
lar mass and star-formation rate; the nature of dust extinction as a function 
of luminosity; the estimate of star-formation efficiency for LIRGs and ULIRGs; 
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the mass-met allicity relation as a function of redshift; the compactness of mas- 
sive, quiescent systems; the high degree of turbulence in the ISM of star-forming 
galaxies; and so on). 

There are some definite limitations in the nature of this review. Indeed, the 
description of physical properties was often still couched in terms of how the 
results applied to a specific galaxy sample, assembled using a specific selection 
technique from among the ones described in Section [2j In order for measurements 
of galaxy properties to have true discriminating power among galaxy formation 
models, it is crucial to use samples that are complete with respect to a well- 
defined property, such as rest-frame optical or near-IR luminosity, stellar mass, 
dynamical mass, or star-formation rate. The majority of the samples described in 
this review slice the underlying high-redshift galaxy population in particular ways 
that add an extra layer of complexity if a comparison with a simulated galaxy 
population is desired. We advocate the design of future surveys with the type of 
physical completeness described above in mind. The NEWFIRM Medium-Band 



Survey (jvan Dokkum et al. 20091 ) represents an important step in this direction, 
but an even deeper survey would be desirable, with both rest-UV and optical 
spectroscopic follow-up. 

Furthermore, we point out that many of the results reported here were (by 
necessity) based on small samples - these include mid-IR spectra of the most lu- 
minous sources, measurements of molecular gas content and dynamics, individual 
emission and absorption-line metallicity measurements at z > 2, high-resolution 
rest-frame optical morphological measurements, and AO-assisted emission-line 
maps. While these results highlight truly compelling questions, they also await 
confirmation from samples an order of magnitude larger for a robust comparison 
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with models. 
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Also, by necessity, some of the most intriguing results regarding gas, dust, and 
structural properties are limited to the luminous extreme of the galaxy popula- 
tion. With the steep luminosity functions described in Section [3. H galaxies with 
L < W^^Lq comprise the bulk (~ 75%) of the luminosity and star-formation 
rate density of the Universe at z ~ 2. Even galaxies with L < 10"*^^ Lq may 
contribute ^ 30 — 35% of the bolometric luminosity density at this redshift 



(|Reddv et al. 20081 ). We must find ways of extending our physical studies to- 



wards fainter luminosities - either by using gravitationally-lensed objects, or the 
next generation of instruments and telescopes. These faint objects are impor- 
tant, and perhaps more analogous to the galaxies playing a crucial role in the 
reionization process at z > 6. 

We close by highlighting two important observational challenges, and looking 
towards the future. First, the direct measurement of gas inflow (accretion) and 
outflow (star-formation and AGN feedback) received only passing reference dur- 
ing the course of this review. The ubiquity of galaxy-scale outflows among z >2 
UV-selecte d galaxies has long been known on the basis of rest-frame UV and opti- 



cal spectra (lAdelberger et al. 2005 



Steidel et al. 20ld . 



20031 . 



Pettini et al. 2001 



Shaolev et al. 2003 



2004 ) . yet obtaining robust constraints on the physical prop- 



erties associat ed with these outflo ws (e.g., mass outflow rates) remains a challenge 



(but see, e.g., 



Steidel et al. 2010l ). At the same time, there is no obvious connec- 



tion between these observ ations and the ex t remely popular theoret ical models 



of cold gas accretion (e.g.. 



Dekel et al. 20091 . 



Keres et al. 2009 



20051 ^. An open 



challenge is to observe a "smoking gun" of gas accretion in high-redshift galax- 



ies. Second, the study of galaxy environments at z > 2 is a field in its infancy. 



Galaxies at z = 2 — 4 83 

While overdensities of LAEs have been identified near luminous radio galaxies 



(|Venemans et al. 20071 ). and a small number of apparent protoclusters have been 
discover ed serendipitously during the course of high-redshift galaxy spectroscopic 



surveys (jSteidel et al. 1998 



20051 ). the study of the dependence of galaxy prop- 
erties on environment at z > 2 remains largely untapped. Environmental studies 
will require extensive spectroscopy of mass-complete samples in both overdense 
and "field" environments, and will have potentially fundamental implications for 
understanding the origin of local galaxy environmental trends. Finally, we look 
forward to the truly exciting insights that will be made possible by upcoming 
facilities, including sensitive multi-object near-IR spectrographs on 8 — 10-meter 
class telescopes, ALMA, JWST, and the extremely large ground-based telescopes 
of the next decade. 
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Figure 1: (From lSteidel et al. 20041 ) Two-color (C/„ — G vs. G — TZ) diagram from 
one of the UV-selected survey fields, demonstrating the UV-selection technique 
described in Section [2.21 The green and yellow shaded regions are the z ~ 3 LBG 
color selection windows, while the cyan and magenta regions are used to select 
galaxies at z '^ 2.0 — 2.5 and z ~ 1.5 — 2.0, respectively. 
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Figure 2: Left: (From 



Franx et al. 20031 ) Js — Kg (Vega) color as a function of 



redshift for several different galaxy spectral types, illustrating the DRG selection 
technique. Solid curves indicate single-age stellar populations with ages of 0.25, 
0.5, and 1 Gyr. J 3 — Kg colors exceed a value of 2.3 at z > 2 as a result of 
either the Balmer or 4000 A break moving into the Js band. Dotted and dashed 
curves indicate models with continuous star formation with ages and reddenings 
of 1 Gyr, E{V - B) = 0.15, and 100 Myr, E{B - V) = 0.5, respectively. The 
dash-dotted curve indicates the color evolution of a single-burst population that 
formed at z = 5, and it also satisfies the color criterion above z = 2. Right: 
(From iDaddi et al. 2004i ) Two-color (z — i^ vs. B — z) diagram for galaxies 
in the GOODS field, illustrating the BzK selection technique. Solid triangles 
represent star-forming galaxies at 2; > 1.4, solid circles represent passive galaxies 
at z > 1.4, empty squares are objects with no measured spectroscopic redshift 
but with photometric redshift Zp/iot > 1-4. X-ray sources are circled, solid squares 
are objects at z < 1.4, and star symbols indicate Galactic objects. The solid and 
dashed lines indicate the regions of BzK color space used to identify star-forming 
and passive galaxies at 1.4 < z < 2.5. 
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Figure 3: (From iGronwall et al. 20071 ) Filter bandpasses used for narrowband 
selection of LAEs at z = 3.1, where a spectrum of a typical z = 3.1 Lya-emitting 
galaxy (overlaid for comparison) would be detected as having a red broadband 
minus narrowband color. This technique preferentially selects objects with bright 
line emission (and, often, faint continua). 
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Figure 4: Multiple determinations of the z ^ 2, ^: ^ 3, and z ^ 4 r est-frame UV 
luminosity functions. Left: z ~ 2: (From 



Reddv fc Steidel 20091 1 Points from 



Reddv &: Steidel (20091 ) are shown as sohd blue circles, while the results from 
other groups a re indicated as in the legend, inclu ding survey areas. Mid dle: 
z ~ 3: (From lEeddv fc Steidel 2009i ;i Points from iE,eddv fc Steidel {2^^ are 
shown as empty squares while the results from other grou ps are as indicated in 
the legend, including survey area . Right: z ~ 4: (From 



Points from 



Bouwens et al. 20071 ) 



Bouwens et al. (20071 ) are indicated as red circles. Results from other 



groups are as indicated in the legend. 
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Figure 5: (From 



Marchesini et al. 20071 ) Rest-frame Optical Luminosity Func- 



tions. Comparison of rest-frame ^-band luminosity functions at z ~ 3. Total 
rest-frame y-band luminosity function of i^-selected galaxies at 2.7 < z < 3.3 
(black dashed line). Also shown are the T^-band luminosity function of "blue" 
(J—K < 2.3) X-sele cted galaxies (blue cir cles), and the z ~ 3 LBG F-band lumi- 



nosity function from 



Shapley et al. (200lh (green stars). Black error bars on the 



blue points represent Poisson uncertainties, whereas gray error bars also include 
field-to-field variations. The gray shaded area represents the la uncertainties of 
the luminosity function of the blue, K-selected galaxies. The inset shows the 
best-fit value and 1, 2, and 3a confidence intervals on a and M* for the blue, 
i^-selected galaxies (blue, solid curve), and LBGs (green, dashed curve). 
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Figure 6: (From ICaputi et al. 20071 ) IR Luminosity Functions. Bolometric IR 
luminosity functions for star-forming galaxies at z ~ 2 (individual detections: 
circles; stacked result for fainter galaxies: square) and z ~ 1 (triangles) in 
the GOODS fields. The local IR luminosity function is shown for comparison 
(crosses). 
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Figure 7: (From ICassata et al. 20081 ) Rest-frame U — B color versus stellar mass 
in six bins of redshift, from z ~ to z ~ 2. The bottom left panel indicates 
the color-magnitude diagram in a protocluster identified at z = 1.61, whereas the 
bottom right panel shows the color-magnitude diagram for all galaxies, regardless 
of redshift. In each panel, the diagonal continuous line indicates a fit to the 
red sequence. The dashed line indicates log(M/M0) = 10.1, the stellar mass 
completeness limit for galaxies on the red sequence. Color-coded symbols indicate 
galaxies in different morphological classes: red, blue, green, and cyan symbols 
respectively represent early-types, spirals, irregulars, and undetected objects. In 
each panel, the inset shows the one-dimensional distribution m.U — B color. A 
clear bimodality in the U — B color distribution is observed to z ~ 2. 
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Figure 8: Star-forma tion History of the Universe. Left: (From 



Hopkins &: Beacom 20061 ) Evolution of the star-formation rate d ensity with red - 



shift. Shown he re is the multiwavelength compilation from JHopkins (200 



(gray points) and 



Hopkins &: Beacom (20061 ) (all other symbols), and references 



therein. All data points have been scaled to a common I MF and dust correction 



The hatched region is the FIR star- formation history from 



Le Floc'h et al. (20051 ). 



The solid lines represent best-f itting parametric forms to the data c ompilation 



from 



Hopkins &: Beacom (20061 ). Right: (From 



Bouwens et al. 20ld ) Evolution 



of the star-formation rate density, based on rest-frame UV luminosity functions. 
At each redshift, the luminosity function is integrated down to Mjjv = —18.3 
AB mag, which is 0.08L* at z = 3. The lower set of points (blue region) shows 
the star-formation rate density determination inferred directly from UV light 
with no dust correction applied. The upper set of points (orange region) shows 
the inferred star-formation rate density using dust corrections inferred from the 
measurement of UV continuum slopes. 
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Figure 9: (Prom 



Maraston 2nnRl ) The effects of TP-AGB Stars on Stellar Popula- 



tion Synthesis Models. Shown here is a compa rison of a 1 Gyr, solar metallicity 



stellar population model from 



Maraston (20051 ) both with (solid, thick line) and 



without (solid, thin line) the contributions of TP-A GB stars. Also shown ar e 



other models from the literature, including those from 



Bnizual fc Gharlot i;20n3l ) 



PEGASE, and Starburst99, as indicated in the legend. Differences in the treat- 
ment of the TP-AGB phase among stellar population synthesis codes lead to 
significant discrepancies in the predicted rest-frame near-IR spectra and lumi- 
nosities of simple stellar populations with ages between 0.5 and 2.0 Gyr. 
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Figure 10: (From iMuzzin et al. 20091 ) Sample population synthesis fits to op- 
tical through mid-IR SEDs, and using different stellar population synthesis 
codes. Shown here are well-sampled multiwavelength SEDs including broad- 
band optical and near-IR {UBVRIz' JHK) and 5piteer/IRAC photometry (red 
circles), as well as GNIRS near-IR spectroscopy (cyan points) for a sample of 
sp ectroscopically-confirmed galaxies at z -^ 2. Also shown are a comparison 



of 



Bruzual &: Chariot (2003i ') (black curve), updated Chariot & Bruzual (blue 



curve), and 



Maraston (20051 ) (green curve) models. 
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Figure 11: (From 



Marchesini et al. 20091 ) The results from various surveys on the 



evolution of the global stellar mass density. Mass-density estimates were obtained 
by integrating stellar mass functions over the mass range, 10® < Mstar/^Q < 
10^^. Red symbols re present the total stellar mass densities estimated from 



Marchesini et al. (20091 ) (i.e., "M09"), where shaded boxes do not include the 
systematic uncertainties, and error bars do. Other estir nates of stellar mass den - 
sities come from the literature, with references as in iMarchesini et al. (20091 ). 



Horizontal dotted lines represent 50%, 25%, 10%, 5%, and 2% of the total stellar 
mass density at z = 0.1. 
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Figure 12: (Prom iDaddi et al. 20071 ) The relationship between SFR and stel- 
lar mass for z ~ 2 star-forming galaxies in the GOODS fields. Only galaxies 
with 24/im detections are included. Left: Star- formation rates are derived from 
UV luminosities, corrected for dust extinction. Right: Star- formation rates 
are derived from the sum of mid-IR and UV-uncorrected luminosities. In both 
panels, filled and empty squares are objects with and without spectroscopic 
redshifts, respectively. The large green squares are the result of the average 
SFR — Mstar relation in GOODS-N based on radio stacking of i^T < 20.5 galax- 
ies in three bins of stellar mass. The blue solid line is the functional form: 
SFR = 200(M,t„r/W^^ Ma;)'^-^(Mf^ vr~^). T he cyan solid lines are the z = 1 and 



z = 0.1 correlations from 



Elbaz et al. (20071 ). The cyan dashed line is a predic- 



tion for z = 2 from the Millennium simulation and semi-analytic models. The 
m agenta star indicate s the location of typical SMGs in this diagram. According 



to 



Daddi et al. (2007i ). star- forming galaxies at 2; ~ 2 appear to follow a strong 



correlation between star- formation rate and stellar mass. This trend is off'set 
towards higher star-formation rates at fixed stellar mass, relative to the corre- 
lations observed among star- forming galaxies at lower redshift. The theoretical 
model shown here for 2 ~ 2 galaxies tends to underpredict the active rates of 
star formation for a given stellar mass. 
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Figure 13: Dus t Extinction in High-redshift Galaxies. Left: (From 



Reddv et al. 2006i ) Dust absorption, parameterized by Ffir/Fi^qq vs. rest-frame 



UV spectral slope, /3, for galaxies at 1.5 < z < 2.6. Filled and open sym- 
bols, respectively, indicate galaxies with inferred stellar population ages of > 100 
and < 100 Myr, for UV-selected galaxies (blue), BzK galaxies (purple), and 
DRGs (red). Filled yellow circles indicate SMGs. The large blue pentagon shows 
the results for UV-selected galaxies undetected at 24/im, using 24/xm stacking 
results. The green filled circles re present the results from X-ray stack ing anal- 
ysis. The solid line indicates the iMeurer. Heckman fc Calzetti f 19991) relation 



Siana et al. 20091 ) 



found for local UV-selected starburst galaxies. Right: (From 

Similar quantities, but for local starbursts (filled black circles), and two z ~ 3 

gravitationally-lense d objects, cB58 and The Cosmic Eye (empty red diamonds). 



Also shown are the 



Meurer. Heckman fc Calzetti (19991 ) relation (solid curve 



and the predicted relation between Lpjfj/XLx^imo ^^id /3 for different reddening 
curves (Calzetti, LMC, and SMC). 
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Figure 14: (From 



Pope et al. 2008l l Mid-IR Spectrum of UL IRGs at high rec 



Pope et al. (;2008l ) 



shift. The composite Spitzer/IRS spectrum of 12 SMGs from 

is shown as the sohd black hne. The purple, red, and green solid his- 

togr ams are composite s from the sample of Spitzer 24;um-selected ULIRGs 
I I 

from iSajina et al. (2007l ). with purple indicating the composite of 4 strong- PAH 
sources, green indicating the composite of 14 weak-PAH, power-law-dominated 
sources with significant 9.7-/xm silicate absorption (rg,/ > 1), and red indicating 
the composite of 17 weak-PAH, power-law-dominated sources with weaker 9.7-^m 
silicate absorption (rg.y < 1). The light-blue dashed curve is the mid-IR spectrum 
of the local ULIRG, Arp 220. All curves have been normalized at ~ 7fira. The 
numbers in the legend indicate the number of sources in each composite spec- 
trum. The prominence of PAH emission in the SMG composite indicates that 
star formation is the main energy source powering the mid-IR spectrum. On the 
other hand, the 24/im-selected ULIRGs indicate a range of spectral types, the 
majority of which appear to be dominated by AGN emission. 
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Figure 15: (Adapted from 



Genzel et al. 20101 ) Relationship between FIR lumi- 



nosity and CO(l-O) luminosity in star-forming galaxies at low and high redshift. 
Gray crosses indicate isolated star-forming galaxies at z ~ 0. Magenta crossed 
squares are merging ULIRG galaxies at z ~ 0. Black crosse d circles are z ~ 1 
star- forming galaxies from the AEGIS field dDavis et al. 20071 ) . Filled blue circles 



are z ~ 2 UV-selected "B X" galaxies (ISteidel et al. 2004 ). Filled green triangles 



are z ~ 1.5 BzK galaxies (iDaddi et al. 201Q) . Red empty squares are z ~ 1 — 3.5 



SMGs (SMG references in 



Genzel et al. 2010l ). In cases where upper- J CO tran- 



sitions were observed (i.e., most z > 1 systems), empirically-calibrated correction 
factors were applied to infer the CO(l-O) luminosity plotted on the horizontal 
axis. Left: Lpm plotted vs. Lco,i~Q- Dashed gray and red lines indicate, 
respectively, the fits to the data for non-mergers (i.e. isolated star-forming galax- 
ies at z ~ 0, z ~ 1 star-forming galaxies, and UV-selected and BzK galaxies 
at higher redshift) and mergers (i.e. z ~ merging and interacting galaxies, 
and SMGs at higher redshift). Right: Lfir/Lco,i-o vs. Leo, i-o- The left- 
hand vertical axis, Lfir/Lco,i~g^ can be re-expressed in terms of the specific 
star- formation rate (in units of Gyr^^), shown as the right-hand vertical axis. 
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Figure 16: (Prom lErb et al. 2006al ) The M^tar - Z relation at z ~ 2. Left: 
Relationship between [Nil] /Ha and stellar mass. Shown here are composite 
Keck/NIR SPEC spectra of the Ha and [Nil] region for the 87 galaxies in the 



sample of lErb et al. (2006al ). grouped into six roughly equal bins of increasing 
stellar mass. In each panel, the mean stellar mass in the bin is indicated, and 
the Ha, [Nil], and [SII] lines are marked by dotted lines (left to right, respec- 
tively). As stellar mass increases, the ratio between [Nil] and Ha increases as well. 
Right: The corresponding relationship between 12-|-log(0/H) and stellar mass, 
based on the empirical trend shown in the left-hand panel. Large gray circles 
represent the averages in each z ~ 2 stellar mass bin, wi th metallicity estimate d 



Pettini fc Pagel f;2nn4l ) 



from the observed [Nil] /Ha using the N2 calibration of 
Vertical error bars show the uncertainty in the [Nil] /Ha ratio, while the addi- 
tional error bar in the lower right corner shows the additional uncertainty in the 



N2 calibration itself. 
of 



'he dashed line is the best-fit mass-metallicity relation 



Tremonti et al. (20041 ). shifted downwards by 0.5 6 dex. In order to com pare 



Tremonti et al. (20041 ). the 



with the metallicities of ~ 53, 000 SDSS galaxies in 
N2 calibration was applied to the low-redshift sample, shown as small gray dots. 
The filled triangles indicate the mean metallicity of the SDSS galaxies in the same 
mass bins used for the z ~ 2 sample. While the SDSS sample clearly shows the 
saturation of the N2 indicator, the more reliable, lower-metallicity bins indicate 



that z ~ 2 galaxies are ~ 0.3 dex lower in metalllicity at a given stellar mass. 
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Figure 17: (From lKriek et al. 20091 ) Morphologies and SEDs of massive galaxies. 
Left: U — B color vs. stellar mass, for a massive galaxy sample at z ~ 2.3 
from 



Kriek et al. (20091 ) with rest-frame optical spectroscopy. The actual NIC2 
image of each galaxy is used as a symbol, in order to indicate the trends between 
morphology and stellar population parameters. Color coding reflects the specific 
star-formation rate of the galaxies. Emission-line galaxies are indicated with italic 
ID numbers, and AGNs are indicated additionally with "(A)." Large, irregular 
galaxies reside mainly in the "blue cloud", while compact, quiescent galaxies lie 
on a red sequence. The ellipse represents the average la confidence interval. 
Right: Stacked SEDs for blue (botto m panel) and red ( top panel) galaxies in 



the 2 < z < 3 spectroscopic sample of 



Kriek et al. (2008bl ). based on rest-frame 



UV photometry and rest-frame optical spectra. 
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Figure 18: (From 



Forster Schreiber et al. 20091 ) Velocity fields from the SINS Ha 



sample, obtained using SINFONI on the VLT. Shown here are velocity fields for 
30 of the 62 galaxies in the SINS Ha survey. Color coding is such that blue to 
red colors correspond to the range of blueshifted to redshifted line emission with 
respect to the systemic velocity. The minimum and maximum relative velocities 
are labeled for each galaxy (in km s^ ). All sources are shown on the same angular 
scale. The white bars correspond to 1", or ~ 8 kpc at z = 2. The galaxies are 
approximately sorted from left to right according to whether their kinematics are 
rotation-dominated or dispersion-dominated, and from top to bottom according 
to whether they are disk-like or merger-like as quantified by kinemetry analysis 



(jShapiro et al. 2003 ). Galaxies observed with AO are indicated by the yellow 



dashed lines. 
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Acronyms/Definitions 

1. LBG: "Lyman Break Galaxy." Star-forming z > 3 galaxy selected on the 
basis of its rest-frame UV colors, which are indicative of a Lyman Break - 
i.e. significant absorption at wavelengths below 912 A. 

2. UV selection: More general use of rest-frame UV photometry to select 
star-forming galaxies not only at z > 3 using the Lyman Break, but also at 
1.4 < z < 2.5 using different rest-frame UV color criteria. 

3. DRG: "Distant Red Galaxy." Galaxy selected on the basis of a red observed- 
frame J — K color, which corresponds to a Balmer or 4000 A break, or 
significant dust extinction, for objects at z ~ 2 — 4. The majority of DRGs 
are actively forming stars, but typically with higher M/L and stellar masses 
than UV-selected galaxies and sBzK galaxies. 

4. sBzK: "Star- forming BzK Galaxy." Galaxy selected on the basis of its B — z 
and z — K colors to lie at 1.4 < z < 2.5 and be actively star-forming. The 
sBzK criteria are tuned to find galaxies with Balmer breaks, ongoing star 
formation, and a wide range of dust extinction properties. 

5. pBzK: "Passive BzK Galaxy." Galaxy selected on the basis of \ts B — z and 
z — K colors to lie at 1.4 < z < 2.5 and be devoid of star formation. The 
pBzK criteria are tuned to find galaxies with prominent 4000 A breaks and 
a lack of current star formation. 

6. SMG: "Submillimeter Galaxy." Galaxy selected on the basis of its powerful 
luminosity in the submillimeter range of the spectrum. SCUBA 850 /im 
selection has commonly been used to find such systems. The prolific IR 
luminosity is due to reprocessed emission from dust (powered either by 
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stars, an active nucleus, or both). 

7. LAE: "Lya Emitter." Galaxy selected on the basis of having strong emission 
in the Lya feature at rest-frame 1216 A. As a result of the selection based 
on line strength, LAEs tend to be significantly fainter in the rest-frame UV 
continuum than other star-forming sources at the same redshift. 

8. LIRG: "Luminous Infrared Galaxy." Galaxy characterized by 10^^ Lq < 
LjR < 10^'^ Lq. Such systems, along with ULIRGs, are much more common 
at z ~ 2 than in the local Universe. 

9. ULIRG: "Ultra-luminous Infrared Galaxy." Galaxy characterized by L/^ > 
W^^Lq. Such systems, along with LIRGs, are much more common at z ~ 2 
than in the local Universe. 



